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of granting £1,000 to the University of London to be devoted to the furtherance of 
research and higher work at University College. After consultation between the 
University and College authorities, the Drapers’ Company presented £1,000 to the 
University to assist the statrstical work and higher teaching of the Department of 
Applied Mathematics. It seemed desirable to commemorate this—probably, first 
occasion on which a great City Company has directly endowed higher research 
work in mathematical science—by the issue of a special series of memoirs in the 
preparation of which the Department has been largely assisted by the grant. Such 


is the aim of the present series of “ Drapers’ Company Research Memours.” 


Kook 


This memoir was written vm 1885, but owing to the cost of reproducing 
the diagrams I am only now m a position to issue it in its entirety. Many 
years ago one scientific society offered to publish the mathematical portion 
without the bulk of the plates, and a technical institute was willing to take 
the graphical work without the mathematics. Revising the memoir after eighteen 
years’ interval with little sympathy for my work of that period I still thought 
wt filled a distinct gap im our technical knowledge, which no researches in the 
mterval had adequately covered. I believe that not only its results, but tts 
methods, which are those of the master, de Saint Venant, will yet be found 


of some suggestiveness. 
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On Torsional. Vibrations in Axles and Shafting. 
By Karu Pearson, F.R.S. 


1. Tue object of the following paper is to deal with certain cases of impulsive 
torsion which seem of considerable practical importance, but which so far as I am 
aware have not hitherto been considered theoretically. I deal in particular with 
shafts and axles. The following assumptions are made: 


(a) The proportionality of stress and strain is supposed to hold for the material 
of the shafting. The limits of safe impulsive loading given in the paper are calculated 
on the basis of stress remaining proportional to strain up to that limit. The material 
must be considered then as having been brought to a state of ease both for positive 
and negative slide, and that between the limits of this state of ease the generalised 
Hooke’s Law holds. This limit to the state of ease is properly the practical efficiency 
of the material, and experiments on its determination for a given material subjected 
to a given process of working and afterwards tested to a given load ought to replace 
the old habit of considering absolute strength and using a large factor of safety. 


(b) The dimensions of the cross-section treated as uniform are supposed small 
as compared with the length of the shaft. At the same time by the use 
of de Saint-Venant’s torsion factors it is not necessary to suppose that the cross- 
sections remain plane. Our results also will be approximately true even when 
the dimensions of the cross-sections are not extremely small with regard to the 
length of the shaft. The cross-section may have any form whatever. . 

(c) The material will be supposed to have a uniform slide modulus p across the 
cross-section, otherwise the distribution of elasticity may be aeolotropic. In the case 
of metal cast or drawn, there is little doubt that » would be different for the core and 
the periphery; at the same time the core is not that part which contributes the 
greater portion of the resistance to torsion, while in other cases the epidermis would 
be turned off. For practical purposes no great error is likely to be introduced by 
treating » as uniform for all elements of the cross-section. 

(dq) I suppose the shafting subject only to torsional stress; the conditions for 
the safe load when in addition to torsional load there is a tractive, flexural, or shearing 
load, must be obtained by the methods indicated by de Saint-Venant in the Memoir 
on Torsion, Chapitre xu. (Savants étrangers x1v. pp. 522—558) or in Legons de 
Nawer (pp. 372—500). 

1—2 
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2. Let «x be the distance measured from one end of the prism along its central 
axis up to any cross-section (area w); let p be the angle through which the section 
has been revolved owing to torsion round the central axis, then for a ‘cross-section 
of any size the corresponding moment of torsion is given by 


na ON 
M= pv TET EE tesa cee cree enna (i), 


where v is de Saint-Venant’s torsion factor. 

If & be the swing radius of the cross-section about the central axis and it as 
allowable to apply (i) to the case of vibrations, we obtain for the case of torsional 
vibrations the equation 

,cy dM dp 
oe dé da dat’ 
where ¢ is the time and p the density. To obtain this equation we have to assume 
that the statical torsional couple M will be produced in each element of the shafting. 
We shall write this equation 


ay _ 7,0 
=F = FR (ii), 
2. BP oe? 
where J pes ears a bake 0 (iii), 


or J represents the velocity of a wave of torsional vibration along the given shafting 
(History of Elasticity, Vol. 1. Art. 1629 or Comptes Rendus, T. 28, pp. 69—72). 


3. The value of the torsion-factor v has been ascertained by de Saint-Venant 
for a wide range of cases in the above cited Memoir on Torsion and in the Comptes 
Rendus T. 87 (pp. 893 et seq.). In the Comptes Rendus T. 88, pp. 142—7, he has 
shewn that for practical purposes we may represent it by the empirical formula 


4, Let o be the slide, S the shear at any point of the cross-section. Then 
S=po. de Saint-Venant in the memoirs referred to has calculated the value of o in 
terms of +(=d6/dx) the angle of torsion, for any point of a great variety of 
cross-sections. 

Let S, be the value of S at which o passes the limit to the state of ease, or the 
material for practical purposes ‘fails.’ I shall call this the fail-limit of the material, 
and the point in the cross-section at which it occurs the fad-point. If we know Sand 
M in terms of the angle of torsion, we can find by elimination a relation between 


them, we thus obtain a value for the maximum safe couple M, in terms of S, 
in the form 


y’ being a coefficient depending on the form of the cross-section. de Saint-Venant 
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(Lecons de Navier, p. 366) tabulates the values of y’ for a considerable number of cross- 
sections. Thus for: 


Circular Section y = "2821, 
Square Section y’ = '2082, 
Square Section with curved sides and rounded angles y/ = ‘2001, 
Square Section with curved sides and acute angles sy’ ='1919, 
Equilateral Triangular Section y ='1755, 
Star with four points y =" 1607, 
while he also gives the value of / for elliptic and rectangular sections of which the 


diameters have a ratio varying from 1 to 8. 
For practical purposes we may take for: 


Elliptic Section y = °2821/n. 
‘ : H) 
Rect ] Aisin 
ectangular Section eg Sec e Tn’ 


where 7 is the ratio of the lesser to the greater diameter of the section, 7.e. n<1: see 
the Memoir on Torsion pp. 351—414*. The former result being accurate to four 
decimal places, the latter empirical and only approximate. 

For a hollow shaft of which the thickness =n times the outer radius 


y/ = "2821 ene = 2921 {/2 (=; . vn)} ‘Pourerialt 
For such shafts there is obviously a great gain in strength. 

The fail-points in all these cases are the points nearest to the central axis of the 
shaft. Asa rule the greater portion of a shaft will be of circular or hollow circular 
section, but in some cases the cross-section near the terminals will be given a different 
form for the purpose of a joint or fitting. In such a case the square shape will be 
better than any given above except the ellipse, but by referring to de Saint-Venant’s 
tables we see that an elliptic section for which the ratio of the axes is even so great 
as 1°8, will give a stronger fitting than a square section of equal area in the ratio 
of ‘2103 to ‘2082, while one in which the ratio of the axes is 1°2, will be stronger in 
the ratio of 2575 to °2082. 


5. The values of » the slide modulus will vary with the material out of which 
the shaft is constructed. Where no better data are available ~ may be taken 2 of £, 
the stretch modulus, the latter being more generally known or at least more readily 
ascertainable. 

In order to keep before us some idea of the numerical magnitude of the quantities 
under consideration, I will calculate roughly the value of J the velocity of a wave 
of torsional vibration in a wrought iron shaft of circular cross-section. In this case 


* I do not understand on what ground Cotterill bases the formula given in the Table, p. 361 
of his Applied Mechanics. 
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we may take p about 10,500,000 Ibs. per sq. in. I take this value from my reduction 
of Kuppfer’s experiments for iron, No. 1: see History of Elasticity, Vol. 1. Art. 1394. 
For a circular cross-section v =k’, and for wrought iron we may take the cubic foot 


to weigh 480 lbs. Then we have 


Ha Bi EI 
a put 480 


= 10,090 feet per sec. about™. 


Thus the vibration, or the torsional couple will be transferred along a foot 
of wrought iron shafting in about yodoq of a second. In the particular case of 
a simple wave in the shafting a wave of length d will repeat itself every 4/10,000 secs. 


i d : : 
or if 7 represents its length every 10.000 * 7° For such a case with the various 


terminal and other conditions of practice, \ may be 40 or 50 or more times J, hence, 
roughly speaking, the wave would repeat itself only 5 times a second, or this is the 
number of times the maximum stresses recur in a second. 


6. Case (2). To one end of a shaft of | feet length a couple N is suddenly applied. 
The other end of the shaft is either built in, or else it needs a greater couple than N 
to overcome the frictional resistance applied to a wheel or other mechanical arrange- 
ment at that terminal. It 1s required to consider the stress at any time after the 
couple is impulsively applied. 

Let py=K. Then the couple V we will suppose applied at c=/; at x=0, P=0 
always. A suitable solution of the equation (11) is 


=A, {cos (Jmt)} {sin (ma)} + at 2, 


which gives uM=k%, =N ife=t 
x 
and ml = (27 +1). 


Obviously w= for all values of « when t=0, and yw will equal nothing for all 
values of #, when t=0, if 


7 =00 P N 
0= fA = 0. 
; = m Sin Mx + Rr" 


Thus A,, is easily determined by Fourier’s Theorem, and we have 


Pa 82 (ody, 7 ae 
v=Rle-33 rq ry 087 3 Ort 1) tin 5 (2r-+1) a. 


* For the practical numerical calculations of this paper it suffices to take J= 10,000. 
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The series of trigonometrical products admits of being broken up into two, which 


allow of discontinuous summation, namely if ¢ lies between + 4 we know that 


wh 


2 i. ee Ens 
| =sing — 5, 8in 386 + = sin 5? — = 8in 19+ 53. 


Hence we can easily deduce that the expression on the right-hand side is equal to 


(- Ly (b—r7) if d lies between 27 — 1 and 2r+1 a 


, - & i eh ; , 
Since cos Jt Bile te {sin (w+ Jt) 5+ sin («— Jt) =} we have the following 


2b Ps 
results : 

(a) (Jt+2) <l, w=0; 

(b) (Jt+a)>l, <Il+a, yap (e—14+ J); 
N 

(c) Jt>(I+a), <(3/—-2), p= (22); 

(d) Jt >(3l—2), < (31 +2), v= > (w+ 31-0); 

(e) Jt > (31+), < (5/—2), w=0; 

(f) Jt >(5t—a), < (51 +2), y= Fle 51+ Ji}; 
N 

(9) Jt> (51+), < (71-2), y= Fe (22) 


(h) Jt>(7l—x), <(7l+2), y= {e+71—Jh, 
etc. etc. 


There is no difficulty now in writing down the general value of w for any value 
of t. It depends entirely on the number of times the ‘head of the wave’ has passed 
and repassed on reflection at the built-in end the section at which we require w. 

Results similar to the above (a)—(h) would have been reached even more readily 
had we adopted the method of solution by arbitrary functions for Equation (ii). 

The Figures 1 (a) to 1(c) on Plate I give the values after any time at any 
point of the shafting of the angle ~, the torsional couple M=Kdw/da, and the 
slide-speed dy/dt. Distances along the vertical give the length of any section of 
the shafting from the built-in end, the time as given by Jt is plotted off along 
the horizontal. 
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These figures lead to some interesting conclusions; (a) the displacement from 
section to section is continuous and varies from 0 to (N/K) 2a. (8) The torsional — 
couple is discontinuous and takes only three (all positive) values, namely 0, NV 
and 2N. 

We note that the displacement and the couple will take as maxima-values 
double what they would have done at each point under a couple statically applied. 
This is of course only a special case of Young’s General Theorem. (y) The dis- 
placement speed is discontinuous, and can only take three values, 0, (N/K )J, 
and —(N/K)J, one of which we note is negative. 

Like the torsional couple it has, however, instantaneously only two values, 
one for either element into which the shaft is divided by the alternate 45° sloping 
lines. 

The practical rule as to strength in this Case (i) is evidently the general rule of 
all impulsive stress in elastic material: The impulsive couple must not be greater 
than a moiety of the statical couple which would produce a failure of the material. 
With regard to this rule we note that the maximum couple 2N lasts longest at 
the built-in section ; here it lasts for a time = 2//J, that is while the wave of torsional 
vibration travels backwards and forwards along the shaft. In the case of a wrought 
iron circular shaft of 5 ft. length (see Art. 5) this is about 1/1000 of a second. It 
is possible that the torsional couple lasting only this time would not injure the 
shaft like a statical torsional couple would do when applied for an indefinite 
time. We note that for the next 1/1000th second the built-in end suffers no stress, 
after which time the couple =2N is repeated for 1/1000th second and so on 
alternately. 

It seems to me that not improbably this repetition of the stress, although of 
course to a smaller degree of magnitude as the vibrations subside, will produce 
nearly the same injury to the shaft as if there were a statical load upon it equal 
to 2N for half the length of time. We may note further that if the resistance 
to be overcome lie between Nand 2N, then although the gradually applied (statically 
applied) couple N will not suffice to overcome it, yet the impulsive couple N will 
do so. Examples of this principle are common enough, as in turning a screw into 
a hard piece of wood, or breaking by twisting a stiff piece of wire—a sudden twist 
will often be effectual. 


7, Case (x). After r seconds have elapsed since the application of the impulsive 
couple N, let a second impulsive couple N’ be applied. To draw a diagram of 
the wave of torsional couple. 

All we have to do in this case is to shift Fig. 1 Plate I. horizontally through 
a distance Jz, replace N by N’ and add the quantities in each compartment. 
The result will depend on where Jr lies between 0 and 4l, 4l/J being the time 
the wave takes to pass twice backwards and forwards along the shaft, or the 
whole time of a complete set of stress-changes. 
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The typical cases are the following : 


(a) Jz lies breveen 0 and 2/, or N’ is applied before the wave has returned to 
loaded section of shaft. 


(b) Jzr=2l, or load N’ is applied at the instant the wave has returned to the end 
of shaft. 


(c) Jr lies between 2/ and 4/, or N’ is applied between first and second 
reflections of wave from loaded end of shaft. 


(d) Jzr=A4l, or load N’ is applied on second reflection of wave from loaded end 
of shaft. This case corresponds exactly to an impulsive couple V+ N’ being applied 
after the first complete period. 


The couples in these four cases are exhibited in Plate I. Fig. 1. (a) to (0). 


These diagrams contain some interesting results. We note again that double of the 
total load couple V +N’ is the maximum stress couple in the first and third cases ; 
we see further that this maximum stress couple is not imposed at each section of the 
rod during the vibrations. In other words the danger of failure is not as in the case 
of statical torsion the same for every cross-section. Thus in the case of (a), only 
sections between «=0, and wx=/—J7/2 are subjected to the maximum couple 
2(N+N’), in the case of (c) only sections between «=0, and «=J7/2-—l. In 
addition, in both cases the time during which the stress lasts at any of these sections 
varies directly as their distance from the section x=/~J7/2. Thus the built-in end 
of the shaft sustains the maximum couple for the longest time and ceteris paribus, we 
may expect therefore failure to take place at the built-in section. If the built-in end 
of the shaft be strengthened by a boss or other mass of metal, the section immediately 
beyond this will be that of greatest danger. 


In case (b), the maximum couple is only the greater of the two, 2N or 2’, but we 
note that after the application of the second impulsive couple, the built-in end of the 
shaft is never free from stress. A singular result is obtained by putting N=’, in 
this case after a time r= 2//J the couple remains uniform throughout the shaft and 
equal to its statical magnitude. Referring to Fig. 1. (c) of Plate I., we note that if it 
be superposed on itself after shifting horizontally through a distance 2/, then dys/dt = 0, 
for all values of ¢ after t=2//J. Thus we have the following proposition : 


If an impulsive couple be applied to a length of shafting, and if at the instant the 
wave of vibrations returns to the point of application after reflection at a fixed end, 
a second equal couple be applied, then the whole shaft 1s brought to vibrational rest, 


N 


each section being subjected to the statical torsional strain \v = zeal (a— ny . 


Of case (d), there is little to be noted, except that after the application of NV’ the 
shaft vibrates as it would do under an impulsive couple V+ NV’. 
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8. Case (iii). A shaft with one end free and the other end fixed (couple not great 
enough to overcome resistance) is subjected to an impulsive couple N, which after 
r seconds is removed. To find the stress in the shaft. | 


All we have to do is to make N’=—WN in the diagrams Fig. 11. (a) to (d) of 
Plate I. We obtain the results Fig. m1. (a) to (d) of Plate I., corresponding to 
the four possible cases of the duration of the load. 


In each case the shaded part of the strip marks the interval during which the 
couple NV is applied. 


These diagrams suggest some rather important practical conclusions. 


(i) Unless the load be removed exactly after a time 4//J (when the shaft remains 
unstrained), the whole shaft (b), or certain portions of it (# and c), will be subjected 
alternately to positive and negative couples equal to +2N. In the case of an 
impulsive couple not removed there is only a positive stress couple =2N. Thus the 
range of stress when the load is suddenly applied and suddenly removed is double the 
range when the load is not removed. 


Now certain experiments of Wohler (see his work, Ueber die Festagkeitsversuche) 
seem to shew that the range of stress must never be greater than the maximum 
statical safe-stress. If this law be true we have the following results : 


For statical torsion the load-couple must never be greater than M,. 


When the load-couple is suddenly applied but not suddenly removed, it must 
not be greater than M,/2. 


When the load-couple is suddenly applied and also suddenly removed, it must 
not be greater than M,/4. 


It is possible, however, that the results of the very careful experiments upon 
which this law of the ‘range of stress’ is based have been misinterpreted owing 
to a not very clear understanding of the effect of the time element brought out 
by the diagrams above. I shall return to this point later. 


(ii) In Fig. 11. (a), (0), (c) there are times at which the stress-couple vanishes at 
every point of the shaft with the possible exception of a single point of discontinuity. 
The times at which this singular distribution of stress arises are as follows: 
In case (a) and (c) when t=2l/J+7/2. In case (b) when t=3//J. In order to 
ascertain a little more clearly what takes place, we must return to Fig. 1. (a) to (c) 
of Plate IL, we must replace ¢ by ¢t—7, and N by —N, then shifting the figures 
horizontally through Jz, superposing them on the original diagrams and adding 
the quantities in the same spaces we obtain the following results : 
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Case Time ne eu Seay Shaft Value of y Value of dy/dt 
= : loaded section 0 —2SN/K 
SS ut lial al aa { built-in section 0 0 
(0) 31/J built-in section whole shaft 0 —2JN/K 
; a loaded section 0 - 27 N/K 
(c) ai ee viiy bit NA eee section 0 0 


In case (d), dyp/dt as well as % vanishes for all values of x when t= 4l/./. 


Thus we see that although the torsional angle s and the stress-couple are 
zero, yet the slide velocity —dy/dt through the whole or a portion of the shaft 
is not so, so that immediately stresses are again brought into play. 


9. We are now in a position to deal with the important case of repeated 
loading. 


Case (w). A shaft built in at one end and free at the other is subjected 
to a torsional couple N for r seconds. This couple is again imposed for + seconds 
after 7’ seconds have elapsed. The process of loading and unloading being repeated 
every t+7’ seconds, it is required to find the nature of the stress produced. 


It is well known that a wave of vibration set up in an elastic material does 
not continue to move backwards and forwards across the piece indefinitely. The 
energy of vibration rapidly becomes less and less owing to the viscosity of the 
material, and within a very short time it is practically at rest. Hence the 
diagrams such as Fig. 11. (a) to (d) above are only theoretical representations ; the 
successive maxima of stress-couples get less and less, till within a comparatively 
short time after the removal of the load the maximum stress-couple produced 
by the motion of the wave is a negligible quantity. Let t, be the time it takes 
for a wave of torsional vibration to be dissipated. Then the solution of the problem 
of repeated load will depend upon a determination of ¢,, and its relation to the 
duration of loading 7 and the interval of no loading 7’. If 7’>¢, then repeated 
loading will have no more effect than a single loading, if 7’<¢, then the maximum 
couples produced by different loads may occur at the same section at the same 
instant, and in this case the repeated load may produce a very disastrous effect ; 
the only case in which 7’ whether > or <¢,, will have no importance, is that 
of case (d), when the time of loading += 4l/J; under these circumstances repeated 
loadings, whatever the rapidity of their succession, will have no more effect than 
a single loading. This case in which the period of loading is equal to, or an 
exact multiple of, the wave period, is very unlikely to occur in practice, and 
as it possesses no features of special interest may be dismissed without further 


consideration. 
: a 
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In experiments on repeated loading I do not think that the importance 
of noting the exact relations between 7, 7’ and t, has always been recognised. 
It is obvious that the quantity ¢, is an essential feature in the working strength 
of a material subject to rapid and discontinuous changes of load; such for 
example as must frequently occur in shafting driving machines, which have to 
overcome discontinuous resistance, or possibly in the impulsive couple produced 
in the axles of railway waggons by passing the junction of two lengths of rail. 
I shall return later to such cases of revolving shafting. Experiments to determine 
the absolute value of ¢, would be impossible, but there would be little difficulty 
in determining how many times the wave runs backwards and forwards along 
the beam after the removal of the load before the maximum couples become 
negligible*. Examining Fig. 1. (a) and (b) of Plate I. we see that this is prac- 
tically the same problem as determining the manner in which the amplitudes of 


v( = RK or 21 at the free end) diminish. 


In shafting of fair section such as would be used in machinery the terminal 
conditions are generally such that the distance between maxima of the stress- 
couples which may be taken to correspond to the wave-length is many multiples 
of the length of the shaft, hence a considerable interval elapses between the 
maxima, and viscosity has more opportunity of influencing the maximum than 
if as in the case of Art. (6) the wave-length is only four times the length of the 
shaft. Thus in such cases we shall be within the bounds of safety if we take the 
greatest of the first two or three maxima. It is possible that the effect of viscosity 
is increased when the area of the cross-section is increased, but I know of no experi- 
mental data on the decrease of amplitude of the torsional vibrations of a cylinder 
with a diameter of several inches. In default of any evidence as to the magni- 
tude of ¢t,, we must content ourselves with giving it an arbitrary value and then 
take examples based upon this value for repeated loading. If we do not take 
any account of viscosity it is obvious that any shaft might be destroyed by 
repeating even a very small load-couple at very long intervals a sufficient number 
of times. 


* Professor A. B, W. Kennedy kindly made experiments for me in 1889 on rods of fair size 
(e.g. % diameter and 50” long). The oscillations (and consequently the stresses) due to an impulsive 
couple were of sensible relative magnitude for 18 to 20 secs. after the removal of the load. The 
reduction was, however, such that I think we should have to look for the maximwm maximorum 
among the first two or three maxima. In the case of rotating shafting and axles subject to impulsive - 
stress, the circumstances are such as to make the wave-length a very much greater multiple of the 
shaft-length and hence to increase the effect of viscosity between the maxima. What we need 
are experiments on axles and shafts such as occur in practice rather than on the wires and rods 
of the physical laboratory. Researches on repeated loading require to be undertaken with a full 


appreciation of the importanee of the time interval factor; this can hardly be said of Wohler’s. 
experiments. 
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In Wohler’s experiments of 1870* no data as to the time element are given, 
but as in some cases he gave upwards of 43,000,000 loadings, it is obvious that an 
experiment of this kind would have lasted longer than a yeart, had the loading been 
made only once ina second. But it is not probable that for the purpose of his experi- 
ments he kept his steam-engines running day and night for a length of time like 
this. We may I think assume that Wohler’s loads were repeated in all cases before 
the stress-wave set wp by the previous loading had ceased to be of importance. 
Indeed it seems highly probable that several successive loadings were producing 
effect at one and the same time; hence I am inclined to doubt the generality 
of Wohler’s Laws and to believe they will ultimately be found to be only a 
particular case of more general principles involving the consideration of the time 
element as an important factor. 

I shall turn now to special cases which will indicate the general principles referred 
to above :— 

10. Turning to the diagrams Fig. 11. (a) to (c) of Plate I., we note that the 
maximum stress has its longest duration at the built-in end and has the two values 
2N and —2N of opposite sign ; these may {(a) and (c)}, or may not {(b)} be separated 
by intervals of zero stress. It will be sufficient then to consider the stress at 
the built-in end and investigate when the maximum stress will accumulate owing 
to repeated loadings. We have only to shift the base line of any one of the 
diagrams through a space Jt’=J(7+7’) to the right to get the effect of a second 
loading. 

(a) We will deal first with only a double loading. 

Now in order that an accumulative effect may arise, the first 2N of the first 
loading must fall when shifted through J¢’, on the second 2N of the unshifted 
loading, or we must have 

(a) Jt’ >4l—Jr< 4+, if J7< 2l 
(b) Jt’ >2l < 61, if Jr= 21 
(c) Jt'>dJr < 61, if Jr > 21 <Al. 


Let m be the fraction of a second the wave takes to run in one direction along 
the bar, then m=l/J, and we have the following conditions for accumulative 
effect : 

t’ lies between 4m-+7, if r<or=2m 
es _ + rand 6m, if 7>2m and <4m. 


It is further necessary that ¢’ be <¢,, the time for relative subsidence of stress- 
wave, or this period must always be greater than 4m—r. If for example 7=2m 
then it is necessary for t, to be > 2m or the wave to run without diminished effect 


* Ueber die Festigkeitsversuche mit Hisen und Stahl, Berlin, 1870. 
+ One experiment with 132,000,000 loadings is said to have lasted 3 years, 
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at least once backwards and forwards along the bar. If 7 be small as compared 
with m, it is necessary for the wave to run at least twice backwards and forwards 
for the second loading to accumulate. I imagine that this condition for ¢, will be as 
a rule fulfilled, and that it is the interval ¢’ between the commencement of successive 
loadings which practically determines the question of accumulation in the case of two 
loadings only. 

We may then conclude that : 

If ¢’ satisfies the relation (i) the maximum stress in a twice repeated loading is 
4N, i.e. double that of a single impulsive loading, four times that of a gradually 
applied statical load. 

Were Wohler’s law correct as to range of stress we should have a range from 
4N to —4N and the safe load would be only one-eighth of a statical safe load. This 
I think is very doubtful. | 

(b) Supposing now the load is applied a third time, then these maximum stresses 
will accumulate provided the condition (i) {modified in case (c) to Jt’ >J7 < 81—Jz} 
is satisfied and in addition (ii)......... rece 

This connotes, if 7 be small compared with m, four runnings of the wave 
backwards and forwards along the shaft without its intensity becoming negligible, 
and if t=2m, two similar movements of the wave. The safe stress limit now 
becomes 1/6th of the statical load on our theory and 1/12th on Wohler’s theory. 
This great reduction of safety is nowhere manifested in Wéhler’s experiments. 
Hence we are forced to the conclusion that the times at which Wohler repeated 
his load and unload did not lead to an accumulation of stress of this overwhelming 
character. 

It is of course obvious that except in the case of an extremely long shaft, 
we have assumed durations for the load and unload of such smallness that they 
are hardly likely to occur in practice,—t’ as a rule will be many multiples of m. 

The general condition will be that ¢’ should lie between 4p x m+7, supposing 
the time 7 to be still very small and p any integer. But these special cases are 
sufficient to bring out how the maximum stress is dependent on the durations 
of load and unload and so the incompleteness of Wéhler’s laws which do not contain 
these durations as a necessary element. In my own mind there is little doubt that 
these durations in Wéhler’s experiments were such as to lead to no accumulation of 
stress due to the first and second loadings. Of course it is possible that the first and 
third (or a later load) may have done, but it seems idle to investigate the question 
without some details of the time data. 

Indeed I have not yet found a method for dealing generally with the question of 
accumulating stress for perfectly indefinite values of ¢’ and r. At the same time 
a system of parallel lines on the drawing board enables one readily to find the 
solution for any definite case. Ten or more loads will be found in a great variety of 
cases to produce no accumulation, while other cases are easily constructed in which 
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the stress would accumulate enormously were it not for the viscous property of the 
material. Hence it seems to me that the laws obtained by Wéhler depend upon 
the rate at which he drove his engine, and that two or three repeated loadings rightly 
applied might have destroyed the material where several million with a different 
period completely failed. 


11. (c) Let us now consider the case where alternate loadings have opposite 
signs. We have in this case only to change the sign of N and shift the diagram 
through the interval ¢’. For generality let us write NM=—WN’. The slightest 
inspection of the diagrams corresponding to Art. 8 shows that, if the interval of 
unload be small as compared with that of load, we certainly get accumulation of 
both positive and negative maxima-stresses, whereas the same time intervals for a 
repeated load of the same sense do not give accumulation. Thus suppose the load 
N can be safely applied to the bar, that is 2 is less than the elastic limit, and let 
this load be continually repeated without accumulation so that it still remains safe. 
Now suppose the load WN, to be alternated with a load —WN,. Then the total 
positive stress, if accumulation takes place, will be 2 (V,+N,) which must be less 
than the elastic limit for safety. This, I think, is the key to Wéhler’s law; it is 
deduced from the ordinary principles of elasticity and involves no recondite charac- 
teristic of repeated loading. We see that the ‘range of stress’ is merely another 
word in the cases considered for 4 {maximum positive stress}, 7.e V,+ NV, and N+0; 
and that the ‘accumulation’ of the maximum-stresses in the case of alternate 
directions of loading under the conditions of Wéohler’s experiments is the probable 
explanation of the law propounded by him. ‘The fact that the law can be explained 
by the ordinary principles of elasticity does not of course detract from the importance 
we must attribute to its discovery. As an example of the remarks in this paragraph 
I give Plate II. Fig. 1v. (a) and (b), two diagrams of a load N, followed after a short 
interval by a load NV’. If N and N’ are of the same sign (a) there is no accumula- 
tion, but there is accumulation if they are of opposite sign (bd). 

We note that in the last diagram Fig. tv. (b) we not only get the accumulated 
stresses +2 (V+N’), but also the stresses +(N+2N’) and +(2N+N’), which are 
greater than any which occur with the repeated loads N and N’ of the same sign. 
Similar results involving, however, different conditions for the time intervals may be 
obtained from diagrams Fig. 111. (b) and (c). 


12. (d) As tending to throw still more light on Wéhler’s experiments I give 
diagrams for the maximum stresses at the built-in end of the shaft for four 
alternate loads NV and WN’, and for four alternate loads N and —WN’ [Fig. v. (a) 
and (b) Plate III.}. Underneath these diagrams are those for six alternate loads 
in the same two cases [Figs. (c) and (d) Plate III.]. Notwithstanding that these 
are only special problems, they indicate that very curious results may be 
obtained by paying regard to the number, duration and sign of successive loads. 
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Thus we find 

3 Two loads NV and NV’. Maximum couples = +2N or +2N’. 

(i) hee loads N and —N’. Maximum couples = +2(N+N’). 

(ii) Four loads VN and N’. Maximum couples= +2(N+N’), the same as in 
two loadings of opposite sign. These maximum couples however do not occur 
frequently and have only small duration. It is very possible that one or two 
only will have occurred before viscosity destroys the effect of the first load 
[Fig. («)]. 

Four loads N and —N’. Maximum couples = +2(N+N’) as in the case 
of only two loadings, but these couples are of much greater frequency and 
duration. In addition series of positive and negative couples occur at the built-in 
end, that is the couples do not so frequently change sign as in the case of four 
loadings of the same sign [Fig. (0) |. 

(iii) Six loads N and N’. Between the fifth and fourth load there is a small 
interval in which the maximum couple takes the value 2(N+WN’), but after the 
sixth load the maximum couples are +2N or +2N’. In other words: a load 
repeated six times is less likely to destroy the shaft than of wt had only been 
repeated four times [ Fig. (c)]. 

Thus to increase the number of times the load is repeated does not necessarily 
decrease the safety. It is thus conceivable that four loads judiciously repeated 
might destroy a shaft, while 43,000,000 applied at stated intervals might leave 
it still intact. Because a piece of material has withstood millions of repetitions 
of the same load with a definite period of rotation of the testing machine, it does 
not follow that it could withstand the same load otherwise repeated and applied. 

(iv) Six loads NV and —N’. Maximum couples = +2(N+W’), as in the case 
of only two or four repeated loads, but the periods over which these couples act 
have been much extended and the danger to the shaft is correspondingly increased 
[Fig. (d)]. 

(v) By superposing Figs. (c) and (d) upon themselves with an interval 
=Jt’, between the first load of the shifted and the last load of the original 
figure we can obtain diagrams for 12 repeated loads. This in most instances 
will probably exceed the limit during which the stress-couple due to the first 
load remains appreciable. In neither case does there appear to be any further 
tendency to accumulate stress, but rather to diminish the duration of the maximum 
couples. We may, I think, then fairly conclude that for this particular problem 
2(N+N’) gives the maximum couple in both cases of the same and of alternate 
sign in loads, but it occurs so rarely and for such short duration in the former 
case that not improbably +2N or +2N’ would be a sufficient limit in practice 
to take for the maximum couples. 

(e) The whole of our work and results with regard to repeated loading has 
been deduced from the solution of equation (ii) in Art. (2). But if # be the 
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stretch modulus of the material of a bar and we put J* in that equation equal 
E/p, then, if for w we read z the displacement, and for N read T an impulsive load, 
the whole of our results will/apply to a bar suddenly loaded and then unloaded, 
one end being supposed fixed. Our diagrams thus give results for repeated 
loading of a bar subject to longitudinal tensile or to alternate longitudinal tensile 
and compressive loads; we have only to put for N and N’ the alternate tensions 
T and T”’. 

It will only be necessary to give J a somewhat different value. For iron 
we may take roughly H=27,000,000 lbs. per sq. in., and taking a cubic foot 
as in Art. (5) to weigh 480 lbs. we have: 


J=16,150 ft. per sec. about. 


Our remarks then on repeated loading apply as much to Wohler’s tensile 
as to his torsional experiments, for both of which, as for flexure, he states the 
same general law. 


13. Although it has not been possible to deduce a general law regulating 
repeated loading for all cases we are still, I think, in a position to criticise 
Wohler’s statements, provided we may assume the phenomena observed by him 
to be due solely to the already well known properties of elastic material. On 
p- 6 of his work he writes: 

Das von mir endeckte Gesetz, dessen allgemeine Giiltigkeit fiir Eisen und Stahl durch diese 
Versuchs-Resultate gepriift werden soll, lautet : 

Der Bruch des Materials lisst sich auch durch vielfach wiederholte Schwingungen, von denen 
keine die absolute Bruchgrenze erreicht herbeifiihren. P 

This is perfectly true, provided the “ Bruchgrenze” here denoted is that 
which would destroy the elasticity of the material by a single sudden loading. 
A single sudden loading produces of course double the stress due to a gradual 
statical loading. The rupture owing to vibrations is due to the effect of stress 
in certain cases of repeated loading tending to accumulate. 


Die Differenzen der Spannungen welche die Schwingungen eingrenzen, sind dabei fir die 
Zerstérung des Zusammenhanges massgebend. 

Die absolute Grésse der Grenz-Spannungen ist nur insoweit von Einfluss, als mit wachsender 
Spannung die Differenzen, welche den Bruch herbeifiihren, sich verringern. 


I believe these last statements to be only partially correct. The first sentence 
I think expresses something near the truth in those cases of alternate loadings 
with opposite sign, in which +2(N+WN’) or +2(7+7") remain the maximum 
loads. Here the algebraic difference is massgebend fiir die Zerstérung des 
Zusammenhangs. But this covers only a certain class of cases in which the alternate 
loadings have opposite signs. It is easy to arrange a series of loadings in which 
the stress shall continue to accumulate above +2(N+N’) or +2(7+7"). It is 
probable that Wéhler’s machine was not, or could not be, run at such a speed 


as to produce any such series. 
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Further, I do not think this sentence even if qualified by the second one, 
expresses the law for alternate loadings of the same sign. In no case that 
I have come across does 2(N’~N) or 2(7’~T) express the maximum load, 
N and N’ being of the same sign. We are thrown back on the case that the 
repeated loads do not produce more effect than the single loads (7.e. no accumulation, 
maximum-couples + 2N or +2N’), or else the maximum stress appears to be the 
same as if the loads had alternate sign {7.e. =2(N+N’) or 2(7+7")} but is 
of very short duration. Of course it is possible also in this case to so arrange the 
loadings as to obtain far larger stress-accumulation. The times suitable for such 
loadings appear to be rare in practice, or we should more frequently hear of 
a small repeated load destroying a structure. I see no reason, theoretically, why 
if 2N approaches the elastic limit of the material it should be necessary to make 
(N —N’) smaller, as seems indicated in the last sentence of Wohler’s law. Wohler’s 
result, however, might be obtained if the effect of his machine instead of giving 
alternate load, N, 0, N’, 0, WN, etc., practically made the loads continuous, J, 
N’, N, ete. In this case, supposing N>N’, we should really be dealing with 
a continuous statical load N’ and a series of loadings N-—WN’, 0, N-N’, ete. 
Without accumulation this would give for the maximum load 


2(N—N’')+N’ 
=2N—N’. 
And as 2N approached the elastic limit it would be necessary to make N’ 
as large as possible, or the difference N—N’ as small as possible, thus it might 


appear that by increasing loading the difference between alternate loads must 
be diminished*. 


* There seems a good deal to confirm this, not only in the form of Wohler’s apparatus, but also 
in the general results of his experiments tabulated on pp. 6 and 7. Wohler it is true is dealing 
with rupture and our theory confines itself to the elastic limit. Still in wrought-iron and steel 
probably any load which fairly exceeds the elastic limit will if repeated often enough “wear out” 
the material, so that I expect for practical purposes Wéohler’s bars were damaged long before the 
conclusion of the often immense number of repeated loads which preceded actual rupture. Let 
T, be the safe elastic limit of tensile-stress, then for a single load 7’ followed by no second load, 
we must have 7'<7'/2; for a single load followed by an equal negative load —7, the maximum 
stress is 2(7’+7') and we must have 7’<T7;/4; for a load 7 followed by 7", (7’<7) we must 
have on the above hypothesis as to the nature of Wohler’s machine (27'—7”")<T7,. Now Wohler 
obtains the following numbers as representing safe limits, z.e. as giving values not equal to destruc- 
tive value of stress; 


7’ = 300, positive load only ; Och pe OCU, 
For Iron iF = 160, alternate positive and negative loads; .. 7,>640, 
7’ =440, 7’ =240, alternate positive loads ; “. Ty>640. 


Hence allowing for effect of viscosity in diminishing previous maximum stress in case of repeated 
loads, we see that these give practically the same limit for 7’. 
T’=480, positive load only; “. T,>960, 
For Cast Steel F = 280, alternate positive and negative loads; .°. 7,>1120, 
7T=800, 7’=350, alternate positive loads ; *. T,>1250, 
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However this may be, I have little doubt, that the best practical measure 
of safety is the ordinary elastic limit, at the same time the ‘accumulated stress’ 
must be taken less than this limit. Each problem in repeated load will require 
to be investigated on its own terms, but so soon as the durations of load and 
unload are known, the problem can be readily solved. Such a rule I think ought 
to replace Wohler’s Law. 


14, Case (v). One terminal of a piece of shafting is fixed, to the other is 
attached at right angles a lever of length b terminating in a mass m. This mass 
by means of a blow (=mV), or otherwise, is started with a velocity V. To 
investigate the stress im the shafting. , 


This case corresponds fairly closely to that of a turn-cock or screw with 
eross-bar at the top, which we attempt to turn by giving a blow at one end 
of the cross-bar. We shall treat the cross-bar as rigid and neglect its inertia. 
There are two sub-cases to be considered. (i) The mass m may be rigidly 
attached to the terminal of the cross-bar. (11) The mass is merely placed in 
contact with it, and leaves it so soon as a tension, not a pressure, is required 
at the point of contact. 

With the previous notation we easily find the following system of equations: 


(i) oe au = for all values of x and ¢ for shaft. 

(ui) mb? ot +K ss when «=0 for all values of ¢. 

(iii) oy V/b, when «=0 and t=e/J, the vanishingly small time which 
the first element of shaft takes to require the spin V/b. 

(iv) w=0 for all values of « when t=0, 

(v) w=0 when «=/ always. 


Now this system of equations is precisely the same as that for the longitudinal 
impact on a bar, which after attempts by de Saint-Venant and by Sebert and 
Hugoniot was finally solved in a finite form by Boussinesq (de Saint-Venant’s 


The divergence here is greater but possibly not greater than would be accounted for by the viscosity 
of axle-steel and the general uncertainty of experiments which stop short of destroying the bar. 


T=500, positive load only; “. T>1000. 

ar Saving Cast Steal | Atvemate positive loadings. 
(Unhardened) 7=700, 1’ =250 *. T,>1150, 
T=800, T’=400 al go= ae 2UU, 
T'=900, 7" =600 “. L,>1200. 


This last set seems to confirm my view of Wohler’s apparatus, 2.e. that it cannot be looked 
upon as suddenly loading, and again withdrawing load in case of positive repeated loads, i.e. 7’— 7” 
must be looked upon as the variable part of the load. 

3—2 
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Clebsch, p. 480k). The analytical results of the latter were converted into readily 
intelligible numerical and graphical forms by de Saint-Venant and Flamant (Comptes 
rendus 1883, pp. 127, 214, 281 and 353). In the following article I make use 
of their calculations interpreted in the language of torsional stress to mark the 
characteristics of the above case. 

We see that in sub-case (ii), the impact will end when dys/dx=0. It is to this 
sub-case (ii) that the curves traced by de Saint-Venant and Flamant can be 
especially applied and it corresponds to a reversal of the couple at the free end 
disconnecting it from the mass m. . 

This seems to me to correspond closely to the nature of a horizontal blow given 
by a body m in a perpendicular direction to a horizontal bar attached to the top of 
a vertical bit of shafting. It would also Transverse section 
correspond to the case where two heavy 
masses A and B, each equal to m/2, 
attached to a cross-bar revolving about 
an axis identical in direction with that 
of the shaft, strike with spin V/b on 
opposite sides this cross-head of the shaft. Indeed this case may be treated as 


covering all impulsive means of putting a shaft in motion im one direction only. 
Interpreting de Saint-Venant’s and Flamant’s results in the language of torsional 
stress we find :— 


(1) That the maximum-couple occurs at the built-in end of the shaft. 
(2) That it occurs before the end of the impact (7.e. before dys/dx = 0, after this 
the motion becomes periodic). 

An examination, however, of the analytical and graphical expressions of the 
Boussinesq solution shows a series of maxima having one maaximum-maximorum 
between the beginning of the impact and the instant when the torsional couple at the 
free end vanishes. When this torsional couple vanishes we find : 

Jt=3:06767 xl; =4°70856 xl; =5:89974x1; =7°41874x1; 
according as y=ls =1/29) =1/4ge= 1/65 
Kl _ palk’ _ second moment of shaft about axis 


where Y= Fm mb? ~ second moment of impulsive load * 


Hence in most practical cases, the head of the impulsive wave will have run three 
and generally many more times along the shaft, before the torsional couple vanishes. 
We may thus fairly assume that if there were any maximum after this which 
equalled or excelled (e.g. with opposite sign) the maximwm-maximorum before this, 
viscosity would in all probability have reduced its magnitude to something less 
or at least to an equality. Thus for both sub-cases we are justified in taking the 
maximum torsional angles and couples before the vanishing of the torsional-couple as 
measurements of the chief strains and stresses due to the impulse. 
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The curves Fig. vit. (a) to (d) of Plate IV. are copied and modified from the 
Comptes rendus (18838, loc. cit.). They give 
(a) The torsional angle % at the free end from the beginning of the motion 
up to the first vanishing of the torsional couple at the free end, for the three values 
of y=1, 1/2, and 1/4. The ordinate gives ¥ =(Jb/V1) ¥. 
(b) The torsional couple M=K dy/dx at the built-in end for the same times, 
and the same values of y. The ordinate gives fA =(Jb/V) M/K. 
(c) The maximum values w,, of % at the free end for values of 1/y ranging 
from 1/2 to 6. The ordinate gives V,,=(Jb/V/) u,,. 
(d) The maximum values M,, of K di/dx at the built-in end for values 
of 1/y from 0 to 25. The ordinate gives f#1,,=(Jb/V) M,,/K. 
From (c) and (d) we see that : 
For values of 1/y > 1, we may safely on 


Maximum value of v=U,= a i 
dy’ 


because the value of Y,, = dee : 2 ay agrees for all practical purposes with its accurately 
3 


calculated value. The latter being represented by the continuous, the former by the 
broken curve. 
For values of 1/y > 5, we may safely take for practical purposes : 


Maximum couple MW, = K NE : 4. if, 


because the value of $#1,,=/1/y+1 agrees for all practical purposes with its accu- 
rately calculated value, the latter being represented by the continuous, the former 
by the broken curve. 


15. Case (vi). One terminal of a piece of shafting carries a wheel whose moment 
of inertia is mB’. A spin w is given suddenly to this wheel (eg. by an impulsive 
couple). The other end of the shafting is ether built-in, or the impulsive stress 
generated is not sufficient to overcome the frictional, or other resistance at that 
terminal. To investigate the stress. 

The system of Equations is now 


(i) oY J a for all values of x and ¢ for shaft. 
(ii) mp = +k as when «=0, for all values of t. 


(iii)  dyb/dt=o, when x=0 and t=e/J the vanishingly small time the first 
element of the shaft takes to acquire the spin w”. 

(iv) Ww=0 for all values of « when ¢=0. 

(v) w=0 when x=/ always. 


* This w must be distinguished from the same symbol used for area of section of shaft, Art. 3, 
and now only indirectly involved in K=pv. This area will be denoted by a for the remainder of this paper. 
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Hence if in the previous case we take b=8, V=fo, then all the conclusions 
drawn for sub-case (1) will apply to this new case, for the equational system will be 
the same for both. 

This case obviously includes the practically important problem of the stress 
produced by suddenly connecting a freely revolving section 
of shafting with a second section, whose frictional resistance 
there is not momentum enough to overcome. 

16. Examples to Cases (v) and (vi). 

Example (1). What is the greatest impulsive velocity 
which can be given to a mass of 10 lbs. connected by a light 
rigid bar; with a short column or shaft, whose other 


i 
wie 
1 
BY) 
» 
i} 
, 


a 


terminal is fixed, without injuring its elasticity, the 
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dimensions being indicated in the figure annexed ? 


The bar is of wrought iron and its mean diameter may 
be taken as one inch. 


{3°6 Cubic inches of Wrought Iron weigh 1 lb.} 


wee ae eae ee 


mx (55) x (P20 x A 


_ second moment of shaft is 3°6 


~ second moment of ball — 12a1G 
= (00037876 = 00038 say. 
We have next to calculate the safe couple for a shaft of these dimensions 


S,=limit to elastic shearing stress = 20000 lbs. say, 4 = 10,500,000, and by equation 
(v) of Art. 4 


M, = ‘2821 x 20000 x wm? x ("25)# inch lbs. 
= 3927 inch lbs. 


Now the value of y above is immensely less than 1/5, therefore we may take for 


the maximum couple Me= aad aie a 1] ; 
bt y 


And since we must have M,, < M,, we have for the limiting value of V 


bM,/ fi “1 
Vis (Vs+1) 


H iK= ceN (5) 1 
ere = 10,500,000 x mw x (5)? x “gz = 1030835 sq. in. Ibs. 
12 x 3927 
Therefore V,=10000 x 04088507 52°3825 ft. per second 


= 8°73 ft. per second. 


Thus the impulsive momentum (10 lbs. x V) suddenly given to the ball must not 
exceed 87°3 speed pounds. 
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The maximum angle sce which the ball is displaced 


SNe 


=5 co radians 


=a! 8-73 J bin Dine: 
= x 10000 x 51°3828 radians 


174 
= a x 51°3828 radians 


001455 x 51°3828 x = degrees 
= 4° 17’ about. 


Example (wv). A propeller shaft is 30 feet long and of 5 inch diameter. It is 
of wrought-iron weighing 480 lbs. per cubic foot. Supposing the propeller to be 
such that its inertia about the axis of the shaft could be represented by a cylindrical 
belt of metal: inner diameter = 30”, outer diameter = 36”, and depth 5”; weighing 
450 lbs. per cubic foot, find the blow which could be given to the propeller 
at a distance of 18” from the axis without injuring the elasticity of the shaft. 

The numbers of this problem are purely hypothetical and may differ very 
considerably from any existing propeller-shaft. The total mass of metal in the 
cylindrical belt would be about 405 lbs. or *18 of a ton. I suppose the great 
majority of propellers would be heavier than this, but it may be doubted whether 
their inertia would be in proportion to the shaft, much greater than that of the 
belt above referred to, which amounts to 111,165 sq. in. lbs. In the next place we 
will suppose the propeller and shafting to be revolving with any given uniform 
spin, one end of the shafting being maintained at this constant spin throughout 
the impulse. In this case to calculate the impulsive stress we may suppose the 
engine end of the shafting to terminate with a fixed cross-section. We should 
then have to add to the impulsive couple thus deduced the value of any couple 
in the shaft due to the resistance overcome by the propeller. 

An approach, however, to the condition of affairs dealt with in the problem 
might be reached if we supposed the rotating propeller to be lifted out of the 
water by the motion of the ship and to receive a blow from the water on 
entering it again. The uniformity of spin at the engine end of the shaft would 
only be maintained by the inertia of the engines and is at best a doubtful 
hypothesis. 

I shall first calculate the statical couple required to destroy the elasticity 
of the shafting in question, we have: 


3 
M,= 8,7 , where a=radius of shaft-section 


= 20,000 x rx ee) inch Ibs. 
= 490,874 inch lbs. 
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The total angle of torsion corresponding to this 
_ 8, length of shaft 
a 5 radius 
_ 2 12x80 
10505 = _2°5 


radians 


ss turns 
~ 3507 
= ‘043654 of a turn, z.e. about 1/25 of a turn 
= 15° 43’ about. 
This is the angle through which the propeller-end of the shaft could be turned 
statically without injuring its elasticity. 
I shall now calculate the value of 1/y. 
{USN _ CSM <5 x 450 


2 
l/y= 


- e x 30 x 12 x 480 


549x99x5x15 
— =! : . 
(2°5)* x 30 x 12 x 16 ails 


This value of 1/y much exceeds 5 and so we may take: 


Kaif 
M,=—$ (+ it, 


n= 9 (a : 


w being the spin impulsively given to the propeller. Putting for M,, the value MV, 

and for K, J and 1/y their values I find »=2°77 revolutions per second, or the 

impulse must be such as to give not greater speed to the propeller than this. 
We find yy, by substituting : 


ED (187117)? 
bn = 79,000 * 2°77 «(eer turns = ‘0349 of a turn 


= 12° 35’ about. 


Thus the angle turned through by the propeller is 20 per cent. less than in the 
case of statical torsion. Let the blow that will just destroy the elasticity be mv, then 
18 x mv=moment of momentum = x moment of inertia of propeller. 

111165 
12 


if mv be measured in foot-second-pound scale. 


18x mv= 


MTT te 


We easily deduce, mv = 8957°2 speed pounds. 
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Now a cubic foot of water weighs 199° lbs, =62°5 Ibs., or 148 cubic feet of water 
falling with a speed of 1 foot per second on a point of the propeller distant 1°5 feet 
from the axis*, would give a blow sufficient to injure the elasticity of the shaft, 
-—always supposing the propeller to absorb the complete momentum of the cubic 
foot of water. 

If the numbers adopted for this problem, as well as its hypothesis of practical 
uniformity of spin at the engine end be even approximately true, the danger to 
the elasticity of propeller shafting owing to the falling of masses of water upon 
the propeller seems very considerable. 

17. Example (wi). A key the head of which is 6” and the shaft 1 foot is made 
of wrought iron, with a uniform cross section of diameter 1” 
in shaft and °5” in the head. 

Supposing an impulsive couple whose moment of momentum 
=1 in foot-lbs.-sec. units to be applied to the head of the key 
and fail to turn it, would this injure the elasticity of the shaft ?— 
We suppose the flexure of the head may be neglected. 


\3 
Let us first find 1/y. The swing radius of head = (5) th 


oe ao ae 
' 


nearly. 
(1/4 x 6x3 at 
ar (°5)* x 12 » LF 


I/y= 


Hence referring to our Fig. vir. (d), we see that M,, = 2°924. 
Thus the maximum couple is given by 


Kae 
M,,=—7 x 2°924. 


4 
Here K=p x 7 where a is the radius of the shaft. 


Now the maximum statical couple is 


TH’ 
M, = S, 2 ’ 
and thus for the limiting spin », we have : 
70 MOE yes 
S; GTX 5 x 7 x 2°924, 
Pa _* x Ly  _ radians per sec. 
eS MeO 2.0 oe 


2 10000 x 12 ih 


pa = 156°3 radians per sec. 
Bids “saws 2924 P 


* Or roughly a cube of water of 2 feet side, falling about 5 feet distance upon the propeller 
would destroy the elasticity of the shaft. Such an impact does not seem improbable in a rough sea. 
4 
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Now wrought iron weighs about 1 lb. to the 3°6 cubic inches. Hence the 
moment of inertia of the head of the key 
_§xmx (d)x6 


T 


in sq. inch lbs. 
ar 3.9 12 8q. Inch Ibs 


Thus we have by the principle of angular momentum 


T 


=1x12x12, 
a 1x 


where © is the spin produced in the key-head by the impulsive couple. It follows 


VAaSeSez 
that O= = we 


radians, or 2 is<@. 

Thus the couple will not destroy the elasticity of the shaft of the key, but 
one very slightly greater (namely, of about 1°066 units) would do so. We see 
that comparatively small impulsive forces are needed to destroy the elasticity 
of such a shaft. 


18. Case (vir). One end of a length | of shafting is maintained at a constant 
spin w [(eg. by being connected with a massive wheel rotating with this spin): 
at the other end of this shafting 1s a wheel of moment of inertia mB which 
together with the shafting is initially rotating with the same spin w. Suddenly 
a frictional couple N is applied to the wheel: to determine the impulsive stress 
in the shafting. 


Measuring x from the end where the friction is to be overcome, we have the 
following system of equations: 


diy _ gt 
TE PoE EESETSSS (i), 
a? d 

ens = K ea, when x=0 for all values of t>0 ......... (ii), 
oY Saw enone tani ii 
yp = When w=l, for all values Of) f . cinesenandh coe (iii), 
dyp 
mes when ¢<0, for all values of w from 0 to l......... (iv) 

v a: p) 

cl 
A te 00100 0.06000 eielee 016 0005.60.01 6100/06 6\u Sie sielsib eee eves sioilela aaininvarainicletanteteninieree (v) 


As a first method of solution, let us employ that of arbitrary functions, 
and take 


p=ot+f(Jt—x)—f(Jt+x—2l) 
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which satisfies (i) and (ii). We have 


d 
EA Ot IS’ (It-2) =f" (Jt+a—2)} 


Tee ee ie Wek ie ee (vii), 
ino —f! (Jt —x) —f" (Jt+a—2l) 
and initially or for t<0 
y= Oator ct wes 
re ys imay take F(t)=0 for Cool eee (viii) 


We have now to satisfy (ii) by substituting (vi) in it. Putting Jt={ and 
K/mB’J*? =p, we have 


£'(O-S" (C21) =p{-F' (QF (¢-2}— Bp, 


FO + RF (O=F" (6-21) —pf (C2) Np. 
Integrating, 
f' ()+2f (0 =P" (6-21) pf (6-21) “PE + ©. 


But when (=Jt=0, /’(¢), f(0), f’(€—21), and J (€—2/) all vanish, hence 
C=0. 
Multiply by e and integrate : 


Fo) = [ ef" (C= 2) ph C21) EF, (pl-1) +0. 
Take the limits between 0 and ¢ and we have 
$ 
f()=er | ef" (L-21) — pf C-2} db— F (pl-1 +e) 


But f(€-2/) and /f’ (€—2l)=0, if €<2I, therefore the lower limit of the 
integral can be changed to 2/. Put y=pl, y=pl=ratio of inertia of shaft to 
that of wheel, and c=n—2y, and we find: 


F(0)= Fe (hn) + f(L=21) = 2° "ef (By) a 


-71 =7-e") + f(C-2)-2€*| “e7(e) do, 


or finally, 


where o=n—-2y, f(c)=f(€—2/), and we have thus the means of determining 


any run of values 2/ of f(¢) from the previous run of 2. 
42 
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19. I have obtained the following series of values of f(¢) by calculation 
from equation (ix). 


(a) f(¢<0)=0, 


F'(C<0)=0 
(6) f(t=9) =x, 0-1-2), 
ies o = ang 1 


(a) F(t=) =F (3) + 3 1) 15 + 4-4) 
+2(n—4y))]x eo}, 


ve — =f" ©) eee +[1+2(n—4y)"]eo™}. 


( 
(¢) f (t= 6) =F (=4)) + ee {ln 6) -7 +17 + 6 (1 6p) + 2 9 6yP 
+3 (9 — by) 
f (t= 5) "(t= 4) + +Rll- —[1 +2 (n —6y) —2 (n — by)’ 
+3 (n= Gy) eat 
f(t=" 5) =f (t= 61) + ie “ {9 —(n —8y) [9 + 8 (7 — 8y) + 4 (n— By)" 
+3 (9 By) eee 
7 (s="3)=0' (61) + +y{-141+4q- 8y)—$ (n — 8y)’ 
+§ (n—8y) Je}. 


20. It is easy to carry the analytical expressions for f(¢) and f’(¢) further, 
but I have not done so on account of the laborious nature of the calculations 
for special numerical values of y, without which the analytical expressions are 
practically useless. In Table I (a) at the end of this memoir I give the values 
of 6(f)=(Ky/N1) f() and ¢'(6)=(K/N)f’ (0) from €=0 to 10/ for the values 
1, 1/2, 1/10, 1/100 and 1/1000 of y. In the construction of this table I have been 
assisted by Newman and Glaisher’s Tables of the Descending Exponential (Camb. 
Phil. Trans. Vol. xu., pp. 145—272), but many exponentials sought had to be 
deduced by additional calculation, as the entries in their tables only go to the 
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second place of decimals in the argument. The labour of determining the values 
of #(¢) and ¢’(¢) much exceeded my anticipations and probably was hardly 
justified by the practical value of the results. At the same time it is interesting 
to note how the curves of spin and stress calculated by means of the exponentials 
give more abrupt changes than those obtained from the superposition of the two 
or three principal harmonics: see Art. 40. The values of the function f() and 
its slope f’(¢) have been calculated and revised with much care, and I believe 
the plotted diagrams will be ‘found sensibly correct for all practical purposes. 
The limitations to the range of values for { only allows in the case of a large 
value of y more than one maximum to be observed, which is not necessarily 
the maximum-maximorum. Allowing for the effect of viscosity, however, it is 
highly probable that the maxima which occur after Jt=10/ do not sensibly 
exceed that which precedes this instant. At the same time definite experiments 
on this point would be of very considerable value. I shall indicate in the sequel 
of this paper, however, how a superior limit to the maximum-maaimorum on 
the hypothesis of no viscosity can be readily obtained. 


21. The graphical results given on Plate V. for y=1/100 and y=1/1000 do 
not differ by visible amounts on the paper from what we obtain by the method of 
a later section of this memoir (p. 57), 7.e. by considering the shaft under statical 
torsion at each instant. Thus ; 

2g’ (£)=1-cos Vy (¢/1+ 1) 
to three and five places of decimals respectively when we compare the tabulated 
values of ¢’(¢) in Table I for y=1/100 and 1/1000 with those calculated for the 
right-hand side of the above relation. The investigation thus shews that for practical 
purposes we may treat shafts for which y=1/100 and upwards, and these are the 
most frequent, by the method of Art. 41. The curves of torsion and stress then 
become such as we have plotted on Plate XI. Fig. xvi. Case (8), or reduce sensibly 
to simple trigonometrical terms. The maximum stress-couple is then =2N, and 
the only noteworthy point is the long interval which elapses before that maximum 
is reached, 2.e. it is not reached until 
Jt/l =1)//y=41°416 for y= 1/100, =99°346 for y = 1/1000. 

Thus, in the latter case the wave must run nearly 100 times along the 
shaft, and during this interval viscosity will generally have had a very sensible 
influence in reducing the amplitude of the wave and so the maximum stress. 

The curves for #(£) and ¢’ (£) are interesting as they exhibit the gradual change 
of form of those functions as they sensibly approach to a simple trigonometrical term, 
when y= 1/1000 and less. 

In the following diagrams and tables I shall only deal with the more theoretically 
interesting if less frequent practical cases of y=1, 1/2 and 1/10. 
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92. The diagrams as they stand exhibit some interesting results, of which I note 


the following : 


(i) At the pulley-end of the shaft, the angle of torsion never changes its sign, 
although it has a greater range of values than at any point nearer the statical end 


(Plate VI.). 


(ii) The angle of torsion has a greater value for y=1 than for y=1/2 
or 1/10 at the 1/4, 1/2 and 3/4 length of shaft, but a less value than y=1/2 at 
the pulley end. Thus to increase the inertia of the shaft relative to that of the 
pulley only lessens the angle of torsion at the pulley-end (Plate VIL). 


(iii) The variations in spin at the pulley-end,—an important point in the 
case of certain weaving machinery—increase as we increase the inertia of the shaft. 
The variation in spin at all points of the shaft is alternately negative and positive, 
and is greatest at the pulley-end of the shaft. 

At the pulley-end of the shaft we have for spin variations : 


30JN _ 307 N : 
y=1, — 865 x —F and +°85 x aye revolutions per min. 
30S N 30SN 
y=1/2, —'65 x By 7 and +°65 x nh ss ” ne 
y=1/10 -°36 x SUES and +°36 x NEES 


ak ah 29 99 ”? 


(iv) The stress-couple has its maximum at the rotated end of the shaft, in 
other words the shaft is more likely to be destroyed at the end at which work is put 
into it, than at the end at which work is taken out,—at the engine end rather than 
at the end thrown into gear. 

At the same time the difference in the stress-couples at the ends of a shaft 
becomes small when the inertia of the pulley is greater as compared with that of the 
shaft: it is still sensible for y=1/10, but hardly for y=1/100. For y=1, however, 
it amounts to as much almost as 2/5 of the applied couple. 


(v) At the pulley-end the stress-couple cannot be negative. It has a very 
slight negative value at a quarter length of span from this end, and from this point 
to the rotated end of the shaft very sensible negative values for large values of y, but 
tends to have only sensible positive values for y=1/10 and less. 


(vi) The dynamical stress-couple largely exceeds the statical couple applied to 
the pulley, and in the case of y being large is considerably greater than its double, 
a value which might possibly have been assumed for it on the basis of similar theorems 
in kinetic stress due to Young, Tredgold, and Poncelet. 

In conclusion I may note that if Wéhler’s views on the influence of range as 
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distinguished from magnitude of stress were correct, then the comparative weakness 
of a shaft under kinetic stress would be still further magnified, owing to the stress- 
couple taking negative values. 


23. Case (vir). One end of a length | of shafting is suddenly started with and 
maintained at a spinw; at the other end of the shafting is a wheel of moment of 
mertia mB*, which together with the shafting is initially at rest, but when revolving 
ut has to overcome a frictional couple N’; to determine the impulsive couple in 


the shafting. 


I shall divide this problem into two parts: (i) to consider what happens before 
the impulsive wave produced by the spin @ given to one end of the shaft reaches the 
wheel at the other, (ii) to discuss the motion after the impulse reaches the wheel. 

This division of the problem is necessary because the couple N’ does not 
commence to act until the impulse has reached the wheel, 7.e. at the time t =//J. 


(i) Since to start with there is only a wave in one direction we choose 


af (Jt—2’) 


where a is measured from the end revolving with spin a. 
Initially aa 0 for all values of x’ from 0 to J, 


=w for all negative values of 2’, 


-, Jf’ (—«a’)=0 unless 2’ is negative or zero, and then it equals o ; 


*, f(¢)=0, unless ¢ is positive, and then it equals ia 


Thus we see that y=0 unless Jt>a’ and then it equals wae t—<’). 


Hence when Jt=/ we have 


oY =a, ow — 5 for all values of a’ from 0 to J. 

This leads us to an important consideration with regard to shafting, namely, if 
one end of it be suddenly started with a spin @, then for safety, it is necessary 
_ that it should be capable of withstanding a greater maximum couple than M=Ko/J. 

This result is quite independent of the length of the shafting and arises solely 
from the impulse given to the shaft by starting one end with the impulsive spin o, 

The importance of a condition of this kind for the safety of the shafting of 


machinery suddenly thrown into gear can hardly be over estimated. 


24. As an example let us calculate the maximum spin which can be suddenly 
given to the end of a shaft 2°5” diameter without injuring its elasticity ; the material 
of the shaft being wrought iron. 
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By Equation (v) of Art. 4, the maximum value of the torsional couple for 
a circular sectional shaft is M,='2821S, /a’, a being the area. 

Further, supposing we may assume the safe tensile load for wrought iron to be 
T, = 25000 Ibs. per sq. in. we should have on the uniconstant hypothesis : 


S =; T= 20000 Ibs. per sq. in. 


Next K=py=pxax4 =), and taking ~= 10,500,000 lbs. per sq. in. as before, 


we have 
2°5 Baas 
10,500, 000 xa x 25) Sie ‘282.1 x 20000 x aVa-—>- - 
8 J 2 
_ 10000 x 12 x 4 x ‘2821 x 20/m 
Hence : 
SES OROO x = revolutions per sec. 


since J =10000 ft. per sec. Therefore # =29°12 revolutions per sec. about. I shall 
investigate in Case (x), the effect on this result of gradually getting up speed, 
we shall see that it is very sensibly modified. 


25. (ii) The result of (i) enables us to measure our time from a new epoch 
and to start with initial conditions, 
eee 
dt” da =F ‘ 
(measuring « from the wheel end of the shaft now), at the instant the couple 
N’ is first applied. Our system of equations now becomes : 


ae 1a » +s g.0 sie'0.6 «es pit eingle Aig ove'rg.oincs alee meg ee acs a en (i), 
mp St = KEN’, when w=0 for all values of > 0 ...ssssse (ii), 
Ww when #=/ for all/ values of ¢.2-2.,,.+--:seeneee eee (iii), 
WY =u, when ¢<0 for all values of a from 0 to /......... (iv), 
WY wd sbeecsssoneceotosbasenSewresendutheptecieey tense t= (v), 


(i) and (iii) give us as solution 


yao (14 25") 4 f(Tt—a)—f(Jt+e— 2) cSt tar eeemeeee tee (vi), 
with the conditions from (iv) and (v) that 
F(C<0)=0 and f’ (C<0)= Oise eee (vii). 
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Substitute (vi) in (ii) and we find: 


mBT LL" ()—F" (0-2) = KY —f'()-F"(E- 2) + GN" 


or, I () + pf () =f" (C— 21) — pf" (—2l) +p 5 a x) 


Integrating and remembering (vii) we find 


F(O+ Pf =F (6-2) = pf (6-20 - (FF) vb 
This agrees exactly with the equation of the previous case (p. 27) if we write 
N’—Ko/J for N- With this modification therefore the series of values for f(£), 
(a) to (e) apply to this case, and the corresponding values of the stress-couple 
will be obtained from the previous results, if we put N’—Ko/J for N in the 
values obtained and add Kw/J to them. 


26. Case (ix). An aale of length 21 terminating in a pair of wheels, each 
of moment of inertia mB’, forms a system rotating with spin w. Suddenly a retarding 
couple N is applied to each wheel: it is required to determine the stress in 
the aale. 


If the shearing stress due to a superincumbent load be added to the torsional 
stress at each terminal of the axle we have a very fair theoretical representation 
of the stress in a railway carriage axle, when the brakes are being applied to 
the wheels, but a still more important application is to plate-rolling mills. 

The system of equations now takes the following form: 


ay ay 
dt? = aL aie. Bie aS Sic aN Wiehe Wala Hele AR in Mie agiaainle a'46)a.e'8 aipe'e 6.6 Simp SU k TMD LAIR W X's 6 aw (i), 
aw dy . ~ 

B de = det both terminals... 006..c-seimers tepecsnsace costs (ii), 

ds : ne 

dn? for all values of a, if t be less than 0 ............... (iii), 

db | 

ae Pare Pim NaNS eh svat. ciel clavalaie wid cc alesele Sls alele Ai wpe ece’m ae Ges eee. 10m, weib.ew GSiAiee he .sieleisld wise (iv) 


Now it is obvious that everything must be symmetrical about the centre 
of the axle, or if we were to take that point as origin of a, the solution 
of (i) must be of the form 


p=ot+f(Jt—a)+f(Jt+2). 
Hence, if we take one terminal as origin, we have 
Watt f(Tt+a—l)t f(TtAI Am) crccesccscrecnscneneenes (v), 


as a proper form of solution. 
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Applying the conditions (iii) and (iv) to (v) we find 
f' (wl) =f" (l-2)=0, 
hence we deduce as before 


f= See f (t= eel ewe Pe (vi). 


Substituting (v) in (ii) we have for #=0 
fh (Se+l)+f" (Jt =p lf (Jt) -f" e+} - FP, 

or integrating and putting Jt={, 
f(C+) + pf C42) =pf (-)—F (C-0 ABE, 


the constant vanishing owing to the relations (vi) when ¢=0. 
Integrating a second time we obtain, if y=p@ as before, 


F(E+)) =7, (1 ewer - (C=) +2] "ef(n-7) sie it era (vii). 


Since f(n—y)=0 for all values from 0 to 2y, so that the lower limit of the 
integral may be taken as 2y instead of 0, (vii) gives us the law of succession 
of the various 2/ ranges of the function f(€) and we note that they closely 
resemble those of case (vu). I tabulate them below. 


=f (= 7) - Fe lln—89) 84 [8 +2 (sy) eo}, 


5 
3 

dh, (< : 5) Sth (< = Bt — F(L-[ $2 (Bye, 

(2) f(b =f) =F (C= 31) + gy O—(9—5y) [9 +4 (7 5y) 


ed Wat d Mc 
i) i (< = i) a Ae 1+[1+2(n—5y)*]eo7}. 
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91 7l\ Nt 
(e) 7(t= mast x i TD Ky WY) TT 8 ST) FZ pete) 
: % Fogo ry) lee Th. 
ats = 34 nah (< = 4 a ate AC Bea) et AC Bea 
+4(n-7y) Je" 7}. 
111 bie ONE 
(sf) f(é = me i Be rae) 97) | Se CR eyya A gay) 
| +3 (n— 97, 
; Hives 9l\ N 
x (c= wee (c= 27) + Ge{-1 411 +4 (n—9y)'- 8 (y—97) 
+8 (y—My)F] oo}, 
Ee laey. 3 WN NL, = _ 
(9) F(o= 1) =F (EH yp) gee ln Hy) 11 +1019) 
+4 (q—1y)+3 (9 —11y)’— 3 (y — Ly) + x (g— Ly) Je}, 
: 131 ; LDN 
if (<= ee (c= op +2 (= My) 4 q- Hy) 
rep pe (Gall y) + te (gy ly) fen ene 
I have carried the formulae one step further than in Case (vii), but they do 


not represent even thus as wide a range of values, as the axle is 2/, and the 


shaft J feet long. 


28. Table I(b) at the end of this memoir gives the values of 


#(Q=F%A(O and o'(O=FP'(0, 


for y=1, 1/2, 1/10, 1/100 and 1/1000, for values of the variable €= —/ to 111. 
In Plate VIII. Fig. x11. (a) are plotted the values of ¢(¢) for y=1, 4, zy, xdv> 
and in Fig. xu. (b) of ¢’(¢), for y=1, 4, x; the others are too close to the 
a-axis to be represented on our scale. Indeed for axles whose inertia relative 
to the wheels is so small, the stress produced by sudden ‘braking’ is also very 
small, z.e. the torsional couple is a very small fraction of the brake-couple. 

Before considering the results obtained, I will first write down the solution 
of the problem on the hypothesis that the axle and wheels form a rigid body. 
We have for the equation of motion 


(mBi-+ yp) GY = —N. 


H a 
ny Sen eat cee 
ence w—w mp (1+y)” 
N re 


¥— Ot TB (+9) "2" 
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el 
But ym? 
Nel ey eee 
Voor hh Wal ee 
Nl y ; 
me cygesa 


Thus the spin is a linear function, and the angle of rotation a parabolic 
function of the time. We should expect therefore as y gets small that the 
values found for W on the supposition of the elastic nature of the axle would 
approach these functions. In our solution we have for the elastic axle 
J —o= 3? [$l (t—a) +4! (Ut+a)} 

oe eee Mee 
— aot = = é 
‘ Koy Y 

For x=0, (or for Jt very great compared with x supposing ¢ and ¢’ to be 

non-periodic,) the terms within the brackets become 


2¢’ (Jt) and 7) 


Hence we should expect for large values of y, that 


2 (Tt) ey. (z) 
Y 4(1+y) Le 
, ies Dé dt 
Gl Fe 2( + ye 
finall PAG) ge eel 
scene ; 4(1+y) P’ 
de SY Salt 
A 2(l+y) 
On Plate VIII. Fig. x11. (a) the values of $ (£)/y are plotted for y=1, 1/2, 1/10, 
and 1/100. In addition the parabola é= Sa is drawn for y=1. It 


coincides extremely closely with the curve € =¢(f)/y, except for values of ¢// less 
than 2°5. The coincidence would be still more marked for smaller values of y: 
We may then say that the angle of rotation is still given notwithstanding the 
elasticity of the axle by a parabolic formula. We can in fact write 


yaor—TT. my (4) + (Are, 


and obtain for all values of t > 2°5 + v/J results giving the angle of torsion sufficiently 


accurately for all practical purposes. For lesser values of t we must use the 
numerical results of Table I (b). 
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Turning to the second figure of Plate VIIL, we see that €=¢'(¢) differs con- 


siderably from €= Sete (the dotted lines). It is true that the curve system 


C+y) 7 
has the same general direction as the line system, but it varies far more than in the 
case of ¢(£) and the corresponding parabolic system. The chief variation is indeed 
the principal harmonic, but the others have for large values of y some influence. 
It would be easy to deduce a formula giving the effect of the principal harmonic 
in addition to the linear term, but the Tables and Diagrams will be sufficient 
to indicate the general nature and practical magnitude of the strain in the axle. 
As we have in the case of the stress-couple to take the difference of two functions 
like ¢’ (€), we obtain no approximation to the stress-couple changes from the linear 
approximation. 


29. It will be found, however, that an empirical formula for the maximum-stress 
couple, which occurs at the wheel end of the axle is,—so far as we are justified in 
deducing one from the small range of values of y considered— 

Max. couple= NV. ieee ; 

Although the stress-couples plotted on Plate X. include only two or three 
maxima, it is still highly probable that the maximum-maxvmorum is owing to 
viscosity included in these values. 


TABLE OF MAxIMA-CoUPLES IN AXLE. WHEEL ENDS. 


y Empirical Formula lst Maximum 2nd Maximum 3rd Maximum 
l NV 864666 V 881732 V ‘784856 V 
(Jt = 20) (Jt = 4°51) (Jt = 8251) 
p | ower | apn | ay | anus y 
5 | SESSA ara oars eal eter 
T00 019802 W dhe 20) Maxima practically equal to these 
Pa 001998 V ae ae occur periodically at intervals of 47/7 


We see that the formula M=N.2y/1+y gives for all practical purposes, when 
y= > 1/10 (as most frequently occurs), the value of the stress-couple. We are now 
able to state a formula for the safe dimensions of our axle, the wheels of which are 
liable to have the brake suddenly applied. 
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Let 2P be the load borne by the axle (as usual beyond its terminals), let NV 
be the couple due to the brake on each wheel. Then the mean shear* due to 
the load, if a be the area of the cross-section of the axle, =P/a, and the shear 
ay 1 
l+y ya 
of equation (v), Art. 4. Hence if S, be the safe shearing stress of the material, 


due to the sudden application of the brake =. , y being the coefficient 
Qa 


we must have: 


P+ NS wee < aS,. 
Ja (1+y)/ 

This formula, as well as the tabular results and the diagrams, apply also to the 
important practical cases of a propeller shaft or run of shafting terminating in a 
rotating mass when the motive power is suddenly withdrawn and the system allowed 
to come to rest by doing work against a constant-couple applied to the rotating mass. 
In this case P will generally be due to the weight of the rotating mass (e.g. propeller, 
pulley, wheel, etc.). Accordingly in the following remarks ‘axle’ may be generalised 
into ‘shaft’ and ‘brake’ into resistance-couple. 


30. With regard to the general results of Table I (b) and Plates VIII. to X. we 
may note: 

(i) That on the application of the brake the uniform retardation of the spin is 
slightly varied owing to the elasticity of the axle. The maximum variation occurs, 
as we should naturally expect, at the centre of the axle. See Plate IX. Fig. x11. (a). 

(ii) The stress produced by sudden braking is very sensible and may for y=1 
nearly reach the value produced by applying the brake-couple in statical torsion to 


the axle. See Plate X. Fig. xtv. (d). 


(ii) For points between the centre of the axle and 3/8 the length of the axle 
on either side (7.e. the mid three-quarters), the stress in the axle is alternately 
negative and positive. The positive values of the stress-couple are very sensible at 


* The result given here as in the numerical examples which follow takes account solely of shearing 
stress ; to obtain the total stress the traction due to the flexure of the axle would have to be combined with 
this shear. If « be the distance from the end of the axle to the point of the journal to which load 
P may be supposed to be applied, the axle, radius r, will be bent by a couple Pe, and the maximum 


traction will be 


Pe 
T= ae RE 
This occurs at all points on the topmost fibre of the axle throughout its length. If 


we have to combine 7 and §. This can be done by de Saint-Venant’s radical formula; or very 
approximately the maximum traction Z” (no longer longitudinal), for an isotropic uniconstant material 
will be, neglecting shear due to P; 
1 
T" = 5— {8Pe+ BN Pe + $8077}, 


(See de St Venant’s Legons de Navier, p. 416.) 


ON TORSIONAL VIBRATIONS IN AXLES AND SHAFTING 39 


the centre of the axle, and if we accepted Wohler’s conclusions as to alternating 
stress wearing out material, it might perhaps be necessary for certain axles to 
consider fully the stresses at points near the centre. 

(iv) It is interesting to note that the maxima and minima follow each other 
with much greater rapidity in the case of an axle to which the brakes are suddenly 
applied, than in the case of a shaft carrying a pulley thrown suddenly into gear. 


31. Applications*. 


I. Each wheel of a truck has to bear a load of 5000 lbs., the axles are made of 
wrought won, the safe shearing stress of which is taken at 20000 lbs. per sq. in. 
Supposing the diameter of each wheel to be 36” and the brakes applied at its ciorcum- 
ference, while the axle has a cross-section of 4 sq. inches; what rs the limit to the safe 
_ ‘braking’ couple which can be applied ? 


Instead of taking P=5000 lbs. in the formula of Art. 29, it seems more rational, 
since the stress is kinetzc, to take it according to Young’s Theorem = 10000 lbs. 
We have then as a necessary condition 


2 
10000 +». 2” — < 0000, 
2°(1+y)y 


N < ‘2821 x 70000 x (1+ 1/y),, N< 19747 (1 +1/y) inch lbs. of moment. 


Or force at circumference < 1097 (1 + 1/y) lbs. 
In most practical cases y will be much less than 1/10 or we must have 


Force < 12067 lbs. 


I have chosen here for the elements values very unfavourable for the safety of the 
axle, namely a very small cross-section and a practically large value of y. 

If we take braking forces of 500 to 2000 lbs., it would appear that we are 
well within the limits of safety. At the same time the sudden application of 
continuous brakes might possibly give rise to much greater values of the ‘ braking’ 
force on the wheels, and in a heavily laden train (7.e. P great) bring the stress 
dangerously near the limit of safety. To investigate this more closely I take the 


following example of a “ theoretical train.” 


Ul. A train is carried on 30 equal wheels, the load is supposed uniformly 
distributed and to be 240 tons; when the speed of the train is a mile per minute the 
continuous brake is applied and brings the train to rest in 200 yards. To find the 
nature of the shearing stress in the aales. 

I will suppose the ‘braking’ force to be 1/2 ton on either side of the wheel, and 
the friction produced between the metals and wheels when the brake is applied to be 


* See footnote on the preceding page. 
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constant up to rest. The ‘principle of work’ gives us at once, if p be the friction at 
each wheel in tons : 


240 x (1760 x 3)? =2 x 80p x 32 x 60 x 60 x 600, 


whence we have 
p= 1°618 tons. 


The total couple applied to each wheel is thus: 


N= (5 ++ st 1-613) x 2 inch tons 


= 32°66 inch tons, 
if 25” be the diameter of the wheels. 

We have next to calculate the value of y and I shall suppose the axle and wheels 
to have the following dimensions. Axle length 4 feet, diameter 4”; the wheels to be 
treated as belts of metal, 25” outer diameter, 24” inner diameter, 4” broad; both axle 
and wheels to be of same material. ‘Then the ratio y of the moments of inertia of 
half axle and one wheel is given by: 

2x12xarx4x4/2 


Ai ee za 
a 4 8 


='0261. 


We are now able to calculate the value of i : eae 


meee ee = 1°64 tons about. 


~ Tax 4 10261 x -2821 


Thus the maximum shearing stress produced by the sudden application of the 
brakes =~ = —_ tons per sq. inch—when we consider torsion only. This occurs 
anywhere on the circumference of the axle at wheel end. 


But there are further stresses which ought to be considered : 


(i) a mean vertical uniform shearing stress due to the weight of the train 
= 8/4 tons per sq. in. 
(ii) a mean horizontal uniform shearing stress due to the friction between the 
metals and the wheels=1°613/4z tons per sq. inch. 
These give a resultant stress =8'16/4a7 nearly in a direction making a slope of 
°202 with the vertical. 


Because the axle is in rotation, it seems to me proper to double this shearing stress 
for its kinematic value, or to take vt: 


16°32 . 
=—7, tons per sq. inch, 


Thus the total shearing stress amounts to 17°96/4a tons per sq. in. or say 
18/47 tons per sq. in. 
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The stress due to the sudden application of the brakes thus consists of two parts, 
(i) That due to the sudden frictional resistance of the metals, 


(1) That due to the sudden torsional couple. The latter in this case forms 
about 9 per cent. of the total shearing stress. Supposing the safe shearing stress 
(elastic limit) to be about 10 tons, we see that in this case there is a considerable 
margin of safety. 


32. Case (x). The end of a length 1 of shafting capable of revolving freely 
ws set gradually in motion, so that at the end of an interval 7, it is revolving 
with speed w, which is maintained after this interval: to investigate the nature 
of the stress. 


My object in dealing with this problem is to mark by a particular example, how 
the interval which a shaft takes to acquire its maximum speed affects the intensity of 
the stress. We have seen in Case (viz) (i), that if the end of a shaft be set suddenly 
in motion with a spin a, it has to sustain a couple= Kw/J which gives a remarkably 
low limit to the spin w. On the other hand, the lowness of this limit need not be 
taken to invalidate our theory, for the end of no shaft can be instantaneously given 
a definite spin; it takes some time, however short, to generate a finite spin in 
all practical cases. That sudden and rapid rotations may injuriously affect a 
shaft is evidenced in certain rolling mills where the motion is very rapidly generated 
and where a special axle piece is put in with the view of its being destroyed and 
rapidly replaced. 

Since our shaft is to rotate freely we have on measuring « from the end set 


in motion: 


WHf (Jt—x) +f (Tt— 24-2). .ccccccsecrconceseeeenees (i). 
Now our conditions for the end set in motion are the following : 
oY = at/r Roe Ta aed DG Pe An Pe ME ORS FACT er paneer (ii) 
SePETOIN TICE ALUCE 652 Ta. coast cere ns nc caovarsequtssatacecapnessaes ass (iii), 
” = 0, when t=0 and for all negative value of : Be na eee ee (iv). 
for «=0 tol 


We have therefore from (iv) 
f(tZ0)=0, 
iv (< zs 0) =0. 
Let us first deal with t <7, then we find from (i) 
ot _ Jf! (Jt) + Jf" (Jt—2l), 
e 
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or it Jt=C 


o£ af (+P (L-2) Bry te eee chr a 


Equation (v) enables us to write down the successive values of /” (2), we have : 


(a) Fo ok 20) 0 
6) f (t= 9)=5-F 
2 


(c) f (t=5)=$- 5 


(d) f (t= )=9-5e- 

OF (Ca)aa3 
endlgeney 

Vins (t= 3 ray) =F — 

oo (eee 5 


This continues until €=J7<2(2r—1)l or <2.2r.1. 


I now propose to write down the values of the stress couple K dys/dz. 


My _ _p (Si—a)+f' (Jt-21+2). 


(a) Jt<zx, fdlh 


dx 
(b) Jt >a <(2l—z), eee 
(c)  -Jt> (21-2) < (21 +2), Ee ey 
(@)  Je> (2t+2)< (ala), GUS EEE 


(e) Jt > (4l—«x) < (4+ 2), a 
(f)  Jt>(4l42)<(6l-2), es 


(9) Jt>(6l—«ax) <(61+2), cp 


dp ow Jt—8l+x 
(h) Jt > (61+) <(8l—2), Isa: = 


We have: 
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(i) — It>(8l-a)<(sl+2), Mo 


? 


dy ow Jt—8l—-ax 
EAPC od Came ROL ee) 

The law of succession is now sufficiently obvious and is represented in 
Fig. vi. (a), Plate Il. We note that (if Jr>2l), the maximum stress-couple 
always takes place at the rotated terminal and has for values eid Fe 

T 

I append, to shew how the shaft assumes the spin o, a diagram of the values of 
dys/dt for Jt<Jzr (Fig. vi. (b), Plate IL). 

We see from these diagrams that although the stresses take a considerable range 
of values, yet that the spin throughout the rod is always of the same sense and 
steadily increases, and that while occasionally the same throughout the rod, it never 


(7) Jt > (81+) <(10/—z), 


differs by more than wl/J7 at any two points at the same instant, and these two 
poimts are then the two terminals or those farthest apart. We easily see that, 
if Jz is considerable as compared with /, the spin becomes practically uniform 
throughout the shaft by the time ¢=7; further, that if Jr be great as compared 
with 2/, the couple becomes practically zero before the rotation is brought up to 
its maximum @. 

On the other hand we observe that if Jr=4rl (r being an integer), whatever 
t be, then the spin will be uniform along the shaft and the couple zero throughout 
when the maximum spin is reached. 

We have now to investigate whether any peculiarity occurs in the stress for 
values of ¢>7, z.e. after the terminal spin becomes uniform and equal to o. 


Now er ae aay 
ai if 


or, it may be considered as the result of superposing an opposite gradual spin 
—w(t—r)/r. Commencing at the end of the interval 7, upon the already existing spin 
wt/r, we have then only to superpose the diagram of the previous article for dy/dx 
shifted through a horizontal distance Jz upon itself, and (# being changed to —@ and 
t to t—7 in the shifted diagram) simple addition gives the required succession of 
stress as in the cases of repeated loading we have previously dealt with. 
Evidently we must have either 
(i) Jr=2r. 21+J7’, where J7’ <2, 


or (ii) Jr=2r+1.21+J7’, where Jr’ <2l. 
In the former case, take Jt=2r. 21+ Jt’, 
in the latter Jt=2r+1.214+ Jt’; 


or we transfer the epoch nearer to the interval t=7. As we are only concerned 
with what takes place after t=7, we shall deal only with that portion of the 


diagram for which Jt=2r.2l or 2r+1.2/ or greater values. 
= 


a ne 
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Case (1). Jr =2r, 214 JS7’. 


Jt’ — 
The stress immediately following 27.2/ now takes the form - + aap or we 


may represent it by the original diagram starting from the origin with Jé’ put 
for Jt; the diagram to be superposed is also the original diagram starting from 
the origin, with —@ for » and Jt’/—Jr for Jt; the shift is to be something 
less than 27. Fig. vit. (a), Plate II. gives the resulting system of stress. 
We notice at once that the maximum positive and negative stress-couples are 
: wo J7’ 
given by Si ey 

Now except in the case where the duration of the load is less than 2/J, 


wo 21 : : 
—-— , hence the maximum stress is reached before 


J Jr 
pe 


the maximum spin at the rotated terminal and equals K sire 


If the wave has run a considerable number of times along the shaft, before 
the maximum spin is reached, 7.e. 27= large number, the stresses afterwards are 


these are always less than + K 


less than K 3 = and are therefore very small. The maximum stresses occur at 
the rotated end. 

If 2r=0, or Jt be <2l, then Jr’=J7 and the maximum stress couples now 
become +Ko/J and take place after the instant at which the maximum spin 
is reached. Obviously then the danger is much increased if the maximum spin 
is attained by the rotated end before the head of the wave has returned to that 
end. The practical rule then is that long shafts which are rapidly rotated should 
be set very gradually in motion. 


Case (7). J7=2r+1.214+J7. 
In this case we note, see Fig. vu. (b), that the maximum  stress-couples 
again occur at the rotated end, but that their values are now 
Fikes se ge Deel 
Sd Te ed ter 


Hence the second pair except when J7’ just equals 27 always gives the maximum 


: 4] — Jr’ 
stress-couple, 2.¢e. +K ae ! gol . IfJd7’=0, Jr=2l, the shaft has to withstand the 


(Jr < 2l). 


large couple + K.2w/J which follows immediately on the maximum spin being 
attained by the rotated end. Thus the most dangerous manner in which a spin 
» can be given to a shaft is applying it gradually in an interval =2J/J; in this 
case the couple amounts to double the maximum which would be reached by an 
mstantaneous application of the spin. 


Whenever the period up to maximum rotation =2r+1.2U/J then the maximum 
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; @ 2 5 oe 
couple is K TF opal and we note that it is very much greater than the value 
wo 1 " é 
K 5. obtained for a period =2r. 2U/J. 


Of course if the period be great (ze. 7 great) these couples may be made 
very small and so quite negligible; when, however, the period of getting up spin 
is only a few times 2//J, the period for sound running along the shaft and back 
again, then the stress-couples produced solely by the setting in motion must be 
taken into account, and it will be important to ascertain whether 7 lies between 
2r—1.21/J and 2r.2l/J or between 27. 20/J and 2r+1.2d/J. 


33. We may illustrate the preceding results by the following numerical example. 
What speed can be safely generated in 1°5 seconds in a shaft 10 ft. long and 
2°5” diameter, the material being wrought iron ? 

In this case taking J=10000 ft. per sec. we have Jr=15000 feet =750 x 2. 
This falls under the first case Jr=2r.2/+J7’. Hence the maximum stress-couple 
_ Ko 21. 1 Ko 
ean 700: J & 

or only 1/750 of the stress produced by a sudden throwing into gear. 


Suppose, however, the time of generating the speed had been 1°502 secs., 
then we have Jr=15020 feet=751x2l. This result now falls under Case (ii) 
and we have the maximum stress-couple 

Ko Al iL Ko 
S Sidia STE Rosi? 


that is practically double of the value in the previous case. As it would be 


impossible in practice to determine whether the speed was generated in 1°5 
or 1°502 seconds, it is obvious that the latter value is to be taken for safety. 
The result shews how the maximum stress varies largely with extremely slight 
differences in the time of getting up speed; a slower generation corresponding 
here to a greater stress. 

By Art. (4), the maximum couple a circular shaft will sustain without injury is: 


M, ='2821 x S,/a® and K=pxax 22) 


For wrought iron we may take »=10,500,000 and S,= 20,000, 
; (2:5 @ iS 
. 10,500,000 x a x 3 x 7X B75 

or w=about 10928 revolutions per sec. 
Thus the safe speed is increased from about 58 to 21856 revolutions per sec. 
in a shaft of these dimensions by spreading the generation of speed over an 
interval of 1°5 secs. If the duration of getting up speed were only so small 


= *2821 x 20,000 xax Jtesy. 
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as 1/100 of a second, still a speed of somewhere about 75 revolutions per sec. 
might be safely generated in this time. A very small time thus suffices to allow 
of any speed likely to occur in practice to be safely generated. 

One point deserves to be noticed: For the same length of shafting the safe 
number of rotations per second varies inversely as the diameter, the time of getting 
up speed being the same. Hence a shaft of smaller diameter might withstand 
better than one of larger a sudden start. Of course the destruction will depend 
not only on the effect of a sudden start, but also on the magnitude of the 
resistance it has to overcome. If, however, a shaft liable to sudden throwing 
in and out of gear, can withstand the torsional resistance, then the smaller the 
diameter the greater the safety. It does not seem improbable that there may 
exist practical cases in which a shaft of smaller diameter would have a longer 
life than one of larger, namely cases in which the impulsive action due to sudden 
alteration of speed, or. even its reversal, is a large factor in the production 
of stress, (eg. in the small bits of axle inserted next the rollers in rolling 
mills with the express object that the shaft, which soon gets worn out, shall 
rupture at this point and not elsewhere). 


34, All the solutions we have yet given to the problem of rotating shafting 
are ina certain sense incomplete; we have only calculated the value of the 
stress up to a certain time, and not its general value at any time after the 
impulse. This result is largely due to the nature of the arbitrary functions we 
have used: the differential-finite-difference equation on which their discovery 
depended did not give us-any general solution. This limitation of time is of no 
importance in practical applications where the important point is to determine 
the maximum stress, which, even disregarding the effect of viscosity, is reached 
very shortly after the impulse and within a few passages of the wave backwards 
and forwards along the shaft. In order, however, to point out one or two 
features of the solutions we have been considering, I propose to discuss Case (vit) 
from the standpoint of that method of investigation which uses ‘normal coordinates.’ 


35, Case (vii bis). One end of a length 1 of shafting is maintained at a constant 
spin w; at the other end of this shafting is a wheel of moment of inertia mB, which 
together with the shafting is initially rotating with the same spin o. Suddenly 
a frictional couple N is applied to the wheel: To determine the impulsive stress 
im the shafting. 


As in Art. 17, we have the system of equations 


mp = = —N, when x=0, for all values of t>0............ (ii), 
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d: 

v=o, Wiel Mote Ale Values OF 4) s.0c.ctete tion eroteeess (iii), 
dis : 
as when ¢<0, for all values of x from 0 to l......... (iv), 
dis 

ai =U) Mavala iietwin atatniatatiatelatyie tei aievelerey ciejaie & eiave sca cies aise cleverare aie eeieia eroeratetachemicitinie (v) 


As a second method of solution let us find a trigonometrical series satisfying 
this system. Assume 


N : 
paott+ F(x—l)+x Ace fc REE Ae OOP: oO ean a (vi). 
Thus we have 
BX _ 7.0X 
ae 4 Fe SESE Cy; 
mp OX = KX, MURDER SOR Me rere setae Cant ags see tts ec ty cs fe nten Sees (ii’), 
dx A; 
TR OSwhen iw =7,etor alli values: of €°%.2..0s lek ees (iii’), 
dy/dx = —N/K, when t=0, for all values of a....... re (iv’), 
d 


H=0, when ¢=0, for all values of x 
The value of x satisfying (1’), (ui’) and (v’) is 
x =A, sin n (l—«) cos nJt, 
where we must determine n from (ii’) and A, from (iv’). 
Substituting in (ii’) we have 
— mB J*n sin nl = — K cos nl, 
or nbtaninb= CK) (MPI my) une tewcssceiscd shes endeee (vil), 


where y has the same value as in Art. 14 (p. 20). The values of n/ are thus the 
roots of the transcendental equation (vii). 


Finally to find A, we have from (iv’) 
N/K= 3A, cos n (l—2). 
By means of (vii) we can easily shew, if n, and n, be unequal, that 
1 
cos n, (/—«) cos n, (1— x) dx =0, 
0 


whence it follows that 
Ni 1 
Z| cos n (J — x) de= A, n| cos’n (1 — x) da, 
0 0 


2Nl1 sin nl 


or A,=y- nl (nl + sin nl cos nl) 
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The full solution is thus comprised in 


2Nl sin nl 


Ns : 
b= ot +(e —I) ai roe MERRIE sin n (l—«x) cos nJt, 


where the summation extends to all values of nl which are roots of 
nl tan nl =y 


The equation (ix) gives theoretically the complete analytical solution, but it 
is obviously of very little value for practical purposes. For every value of y we 
should have to calculate the roots of the transcendental equation afresh, and having 
found these roots substitute in the trigonometrical series and calculate its terms for 
each separate value of # and Jt. The labour involved in this process would be very 
great, and without taking a very large number of terms the solution would be 
nothing like so approximate as the calculations from the arbitrary functions. 


36. The only case where I am aware of a similar calculation having been made 
is that of a beam terminally supported and struck at its centre. This was dealt with 
by de Saint-Venant in 1853—4, when a paper was presented to the Société philo- 
mathique followed by a paper in 1857 to the Académie. (The numerical results are 
given on pp. 553—563, with several corrections at the end of the volume, in the 
annotated Clebsch). The curves for the displacement of the beam at successive 
intervals drawn by de Saint-Venant have never been published and will not now be, 
as he himself felt that they had not been calculated for a sufficiently wide range 
of the ratio of the mass of the beam to the mass of the striking body, nor again for 
a sufficient number of points along the beam. They are represented so far as de Saint- 
Venant’s calculations permit in his plaster of Paris model of impact on a beam. 


Bearing in mind the great labour de Saint-Venant devoted to the particular case 
of transverse impulse treated by him, it seemed obvious that for practical purposes 
the numerical calculations must be based on the more readily dealt with arbitrary 
functions. These arbitrary functions express in fact the sums of the two trigono- 
metrical series obtained by splitting up the terms such as sinn(J—«)cosnJt 
into the sum of two sines. Various interesting identities may thus be obtained 
concerning series in sinnw and cosnz, where n is a root of the transcendental 
equation ni tan nl=y, but their consideration would carry us beyond the limits 
of our present topic. 


One particular case may be dealt with by the trigonometrical method as a means 
of verifying our previous results. Let us suppose y small, or the moment of inertia 
of the wheel large as compared with that of the shaft, we may then throw the 
solution (ix) into a simpler form. There are two distinct kinds of solutions in this 
case for the transcendental equation ni tan nl=y. 


ae ee 
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All the roots but one may be expressed in integer powers of y. I find approxi- 


mately for them if 7 be any integer: 


nl =tr+— i 


ir fr 


( 2). 


On 


y' G 5 


ee 


3 Pr 


which converges rapidly even for not very small values of y. The remaining root 
n, depends upon powers of Jy and I find approximately 


Elie 


nl=y2- a+ 


360" 5040” BOR CORO CCOnGAGOOHOL Goce ( 


xi), 


It remains to substitute these values in the expression for A,. For purposes of 
comparison, I will calculate the quantities involving A,, which occur in the values of 


the stresses at both terminals of the shaft. 


Hence x=/ 


and x=0 


Further we have 


We have 


ee 


eoeeeerene 


ds sin 7l 
K da NN ~2N= mien bcos nl 208” EO -#) cos n Jt. 
sin nl Jt 
M,=N}i ae Be anak oon) oo nl 7 cc eceveves 
sin 2nl Jt 

M,=N}1 —> Paine loeaal oo nl} Serra 
dis Pe eine sin’nl ‘ante 
CE (20) ~~ RK “nl+sinnicosnl> 1 


37. In the following tables I have calculated (i) the values of the first five 


nls for y=1, 1/2, 1/4, 1/6, 1/10; Table II. 
coefficients in equations (xii), (xiii) and (xiv); Table III. 


values are radians, the second are degrees and minutes. 


Nol 


86035 
49° 17’ 40” 


65330 
37° 25° 53” 
48009 
27° 30! 27” 
39725 
22° 45’ 38” 


31105 
it 49° 19% 


* The values here given have been obtained from those 
approximations. 


TABLE II*. For large values of y. 

Ny Nol Nal nl 
3°42562 6°43731 952933 12°64529 
196° 16’ 25” 368° 49’ 50” 545° 59’ 26” 724° 31’ 17” 
3°29231 6°36162 9°47749 12-60601 
188° 38’ 8” 364° 29’ 38” 543° 1’ 12” 722° 16’ 17” 
3°21910 632271 9°45122 12°58623 
184° 26’ 27” 362° 15’ 52” 541° 30’ 54” 721° 8' 17" 
3:19373 6°30959 9°44243 12°57962 
182° 59’ 15” 361° 30’ 47” 541° 0’ 40” 720° 45’ 32” 
3°17310 6°29906 9°43538 12:57432 
181° 48’ 18” 360° 54’ 34” 540° 36’ 26” 720° 27’ 20” 


(ii) the corresponding values of the 
In Table II. the first 


of Equations (x) and (xi) as first 
I believe the radians correct to the last place of decimals, as they were calculated 


7 
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Taste III. For large values of y. 


y Nol nyl Nol Nol nl 
Cy 1:11914 — 15169 04658 — 02167 01239 
Mh eee ‘72986 *14561 04604 02155 01235 
C3 33695 02125 ‘00357 00113 00049 
Cj 1:07014 — 08728 02433 — 01106 00627 
4 <C, 84973 08629 02426 ‘01106 ‘00626 
Ce *32522 00655 ‘00095 00029 ‘00012 
¢, 1:03820 — ‘04697 01242 — 00558 ‘00315 
i <@ 92083 *04684 01241 00558 00315 
C3 23975 ‘00182 ‘00025 ‘00007 00003 
C, 1:02620 — 03213 00833 — 00373 00210 
TCs ‘94629 03207 00833 00373 ‘00210 
Cs 19851 00084 00011 "00003 00001 
C, 1:01609 — 01966 00503 — 00224 ‘00126 
to \% 96734 ‘01965 ‘00503 00224 00126 
Ce "15559 ‘00031 “00004 00001 ‘00000(5) 
rt 4 sin nl 2 sin 2nd 2 sin? nl 
¢* = =——_._~— Ca C3 


~ Onl + sin Qnl’ ~ Onl + sin 2nd’ ~ Onl + sin 2nl 


to 7 places of decimals, the values in degrees, minutes, and seconds were calculated from the values 
of ni to 7 places; a second value differing by 5” was then calculated, both were substituted in the 
equation nitannl=y and a value interpolated between them by a simple proportion to make the 
right-hand side exactly y. In nearly all cases when substituted this value gave the right hand 
a value differing from y only after the fifth place of decimals. The degree values of ni are thus 
in a slight degree more approximate than the radian values. I do not think the maximum error 
can be more than +3”, and this probably only in the column for ,/. 

* The approximate values of ¢,, c,, c, to the third power of y are 


ey = SP fa sy _ (3 0) 


nr? Cr Pr \2 Pr 
2 3y 3 10 
=a {1 a ee — (a = z=) a0 e)e1e/9 0:sieieie.a,einjaiy a6 iho ataleleisieiteriars (xv). 
cz 4y 
a= ga! wat 
I have used these however only for purposes of comparison and verification and calculated the 
c’s from their transcendental values by substituting in them the values of nl given in Table II. 
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For all technical (as distinguished from physical) problems Table III gives 
sufficient coefficients for calculating the values of the maximum stress and of the 
variations of the spin, except perhaps for the case of y=1. If further values 
of the coefficients for this case are required they may be calculated with sufficient 
approximation from the formulae 


é=(—1)4 2 mh ee a) Seles 


Pa vart Cm va’ en 


I add a supplementary table for these additional values: 


Taste IV. High values of nl for y=1. 


y=1 Nl Nel nel Ngl 


C, — ‘00798 00557 — 00411 00315 
C, 00794 00556 ‘00410 00315 
C, 00026 00015 00009(54) 00006 


Any calculations made from Tables III and IV for the maximum stress will 
then have an error of less than 1/2 per cent. 


38. The above large values of y, although theoretically of the greatest 
interest, are not those which most frequently occur in shafting. Such values are 
those lying between ‘05 and ‘005. Accordingly I have added a fourth table, 
Table V, giving a series of additional small values for y and the corresponding 
coefficients for n, and n,/. For practical purposes it is unnecessary to go further, 
and this table enables us to read off very easily the value of the stress at either 
end of the shaft and the variation in the spin, due to the suddenly applied 
frictional resistance. We have for most practical purposes 


M=N 41 — ¢, cos nl —c cos nl = t} 


M=W{1 —C, COS nl —c,' cos nl TH ben Meee whe ete (xvi), 
: Jt 
(Ht). =o- Bs sin n,/ mF 


“%=0 being the terminal bearing the wheel at which the frictional couple NV is 
applied. 
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TaBLeE V*. For small values of y. 


| y Nol a Ce C5 Nyl 

08 | seca age | 101206 | 97376) C1804 ene ie 
06 | ott Sin | 100819 |) -a88bb e elioe yeh et an nome 
03 | o: narogy | 100495 | -99005 “08617 | jan onnaue 
a eee 100332 | -99338 | -o7048 | Jo. 5102, 
015 | gelperdgy | 1:00250 | -99565 | -06109 : ee 
01 | geag ia) | P0016? |. -09607 | a-O4GDe er naa 
008 Say 100183 | -99733 | -04é66 | Ve aeaen 
005 | joards» | 100088 | -s9ss4 | 03533 | 2 IA3le 
003° | grarage | 100080 | "99900 P-oara7 || Veen anne 
002 | ~ a: ggrain | 100088 | 99088" | Rah Aine age 

001 | 4.7502, | 100017 | -99966 | -01581 eye 
0005 | y-r1¢'9” | 100008 | 99983 | -o1118 | , at 

0001 | ge gyo3r | 1:00002 | -99996 | -o0499 one 


* The values of ¢, or c, for n,J—since they are for the above values of y equal in magnitude 
to five places of decimals—have been determined at first by the trigonometrical formula and _atfter- 
wards by the expression ‘202642 x y—:061596 x y*, which was very soon quite sufficient. The values 
of c; for nm, may be obtained if needed from the expression 0322527 — 013072), but for all 
practical purposes their magnitudes are negligible. The values of ¢, and c, for Mf may also be 
readily determined from Equation (xv). But for all technical purposes we may legitimately neglect 
them as well as the values corresponding to n,/ when they are <-03, 

The values of ¢,, ¢, cs for ml have been determined from the trigonometrical expressions by 
substituting in them the value of »,J to 7 places of decimals (i.e. not from the five places tabulated). 


This was continued till the following approximate expressions gave correct results to five places 
of decimals : 


ON TORSIONAL VIBRATIONS IN AXLES AND SHAFTING 53 


39. I propose now to discuss the practical conclusions which may be drawn 


from these formulae, and to plot out the values of M,, M, and (=) for several 
0 


values of y. 

(i) For most practical values of y, ¢, and c/ may be neglected, hence the 
maximum-minimum couples at the wheel which gives out and at the wheel which 
takes in the work are respectively N(1+c,) and N(1+c¢,). Thus we see from 
the tables that the maximum |couple occurs at the end of the shaft which is 
giving out work and is always greater (although often only shghtly) than 24. 
The couple at this end can also change its sign and take a small negative value 
since c¢,> 1. 

On the other hand the maximum couple at the work-receiving end is always 
less than 2 since c, is <1 and thus cannot change its sign. 

(ii). The maximum variations in the spin are +2(JN/K) C,. For many spins 
w and frictional resistances NV, these variations will be small, but in other cases 
it is quite possible that they may be commensurable with w. It is obvious that 
they do not depend upon , but only on the magnitude of N for the same 
shaft and pulley. For a circular shaft of diameter 2a, K=pak? (if a=7a’) and 
the maximum stress-couple has for its limiting value M,=N{1+¢,+¢/+¢,"+ ...}, see 


Art. 42. But M, cannot be greater than S,ak*/a. Thus the greatest possible value 
ae ae Ne fo Rt f c Ronee G| 

—. = s— .—.-.- lu- 

of the variation Ko 3s equal to GesccuT a r to ba ce fw revolu 


In the case of wrought iron when we put »=10,500,000, 


i = To oer lations per second, 


tions per second. 


S,= 20,000, 74=10,000 x 12 in inches, we have 


a being read in inches. 

Hence, if the frictional resistance approaches that which would destroy the 
elasticity of the shaft, it seems possible that, for a large c, (ue. a large y) and 
a small a, the spin at the wheel might even become negative. 

Of course it must be remembered that we have taken a very large value 
for N,—its limiting value. Still it seems important to note the considerable 
variations in spin which can take place owing to a machine being suddenly 


The approximate expressions (x), (xi), (xv) and (xvii) will enable the reader to calculate by the 
aid of the following as many further values of nl and ¢ as he may require. 


: 1 ; if 2 4 + § 

i= in lia 1/i?3? 1/i?r3 L/itr* 1/i®xr6 372 ~ nd | Blin’ ~ a? 
a 3°141593 318310 101321 032252 010266 “001040 004215 002702 
2 6°283186 159155 025330 004031 000642 000016 “001140 000121 
3 9°424779 106103 011258 001194 000127 000001 000372 000017 
4 12566372 079577 006332 000504 | -000040 000000 000162 000004 
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thrown in or out of gear in a shafting. One would expect the effect to be 
distinctly prejudicial in the case of a number of weaving machines on the same 
shaft, where uniformity of spin is so essential to uniformity of work, that one 
powerful engine with a long shaft is frequently replaced by a number of small 
engines each working only a few looms, 


For example, suppose NV to be only 1/8 of the maximum safe couple, and let 
y have the value it has for Case (a) of the following example, then c,=°233 
about and c, may be taken as 1°1, a=2°5/2, then the variations amount 
to as much as 7°7 revolutions a second, and if the shaft were being driven 
at the slow rate of 100 revolutions per minute, the effect of throwing a machine 
into gear might stop and even reverse the spin of the pulley on the shaft. I do 
not know whether any effects of this kind have been observed, the ordinary 
values of y and c, in many cases rendering the variations very much smaller as 
compared with w The effect too, of the usually progressive, if very rapid, appli- 
cation of the frictional resistance would tend to diminish the variation, as we 
have seen in a like case in Art. 33. At the same time it must be remembered 
that the mean value of the spin even for such a short interval as a second is 
practically » and thus positive; the period in such an unfavourable case as (a) 
of the following article for the variation in the spin being only 00674 of a second. 
Thus the negative spin would probably only manifest itself as a scarcely sensible 
vibration of the pulley. 


40. (ui) The point of the above remark will perhaps be brought out more 
clearly by the following numerical example. 

A wrought-iron shaft of 2°5” diameter runs along the side of a shop 50’ long. 
One terminal is started with a speed of 100 revolutions a minute. At the other 
terminal is a cast-iron belt pulley which being thrown into gear is subjected 
to a frictional couple of 20000 inch Ibs. Determine the maximum couples at the 
terminals of the shaft and the variations in the spin of the pulley when it has the 
following dimensions : 


Case (a) Breadth = 4” Case (8) Breadth = 6” 
Thickness of rim = °445” Thickness of rim = 1’ 
Radius = 10” Radius = 20”. 


We shall neglect the spokes of the pulley, and suppose its inertia practically due 
to the weight of the rim. The density of wrought iron to that of cast iron will be 
taken as 480:450 or 16:15=(p/p’). Dealing first with the pulley (a) we have 
for its moment of inertia 


T 


5 X 1664°67 x 4, 


p’ 5 {Lo*— 9555} x 4=p! 
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while the moment of inertia of the shaft 


=p x m(1°25)* x 600 x i) x 5 (1°25)* x 600. 

Hence the ratio y=:2346 nearly. 
In the case of the pulley (8) we have for its moment of inertia 
p! 7/2 {(20)'— (19)} x 6, 
16 ,(1°25)* 
15 ere 33 ue AMET Ts 

It is thus obvious that not quite doubling the size of the belt pulley enormously 
decreases the value of y. 

Before seeking for the values of the constants, I will determine the limit to the 


frictional couple which could be applied statically to a shaft of the above kind. For 
torsion + the couple M is given by: 


and for y the value 


M = pr x moment of inertia of cross-section. 


Hence, if S, be the safe shear limit=yrxa, where a is the radius, we have 
for the maximum couple: 


S,ra* 
M, = oe . 
Taking S,=4/5 7, and T,=25000 lbs. per sq. in. for wrought iron, we have 
(1°25)? 


M, = 20000 x 3°14159 x ous 31415°9 x (1°25)*=about 61000 inch lbs. 


The couple on our pulley is only 20000 inch lbs., but it does not follow that 
we are thus within the safe limits for working stress, because the maximum couple 
applied to the wheel does not measure the maximum couple in the shafting. It 
is necessary that N(1+c¢,) should be less than M,. This is true in our case 
because ¢, is only slightly greater than unity, but it would not be true if the 
frictional couple had been 35000 inch lbs, instead of 20000, a magnitude, however, 
much below the maximum statical couple 61000 inch Ibs, I will now give the values 
of the constants. 


Case (a) Case (8) 
Calculated Interpolated from tables Calculated Interpolated from tables 

C, 1:03604 1:03598 1001461 1:001459 
c, 92544 92552 "997090 997086 
C, 23288 23219 046756 046759 
Cc,” — 04429 — 04423 © — 001773 — 001774 
a 04416 04411 ‘001773 ‘001774 
Cc. ‘00162 00164 0000025 os 
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In Case (a) the interpolated values of the constants have been obtained from 
Table III by taking the difference between 1/4 and 1/6 and using proportional 
parts. In case (8) the difference was taken as the mean of the differences up 
and down from the number (:008) to which y(‘00877) was closest, in order to 
obtain a difference as close to ‘008 as possible and then using proportional parts. 
It will be seen that when y is only sensible in the third place of decimals (= 00877), 
then the interpolated values of the constants agree to five places of decimals with the 
calculated and thus there ought to be no hesitation in using in this case the method 
of interpolation. When y is sensible in the first place of decimals (‘23465) we have 
to work from Table III, where the differences given for y are far greater but the 
maximum error is here only ‘0007 (occurring in ¢,) which we may fairly neglect for 
practical purposes. 


We have also for Case (a), 
nol = "4662673 


for Case (8) 
nl = "0935116 


= 26 3405007 =e) oe 
n,l = 3°21446 nl = 3'1443821 
= 184 01 =180° 930", 
Our numerical results are therefore for Case (a) 
Jt 


M,(=couple at rotated end) = NV {1 — 1:03604 cos n,/ —- + 04429 cos nl} : 


l l 


M, ( =couple at pulley) =e ! — "92544 cos nls! — 04416 cos n,/ ol } 
OQ (=spin of pulley) =o- al {28288 sin 7,1 a +°00162 sin n,/ 7 : 


and for Case (8) 


M,=N 41 — 1:001461 cos nl 54 0017783 cos nit : 


i 


M,=NjI — ‘997090 cos n,/ ] 


001773 cos nl & E 


2. : 
Q=oe- “7 {046756 sin nd + 000025 sin mir} : 


where nJ and n,J have the values given in the table above. The values of 
M,/N, MN and the variation in Q, te. Q—o, are represented graphically in 
Plate XI, Fig. xvi, Cases (a) and (8). In revolutions per minute we have 
2JN/K =1138°32 (assuming J=10000 ft. per sec.) and so in 


Case (a2) Q=100—265°09 sin n,l Jt/l—1°84 sin nl Jt/l, 

Case (8) 9=100— 53°22sinn,J Jt/l —-028 sin nl Jt/l. 
Thus we see a fundamental difference between Cases (a) and (8). In the former 
the spin varies from 367 to —167 revolutions a minute, or the pulley can take 
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a negative spin. In the latter the spin varies from about 153°25 to 46°75 revolutions 
a minute, but is never negative. In both cases we see that very great irregularity is 
produced in the rate of the pulley by throwing it into gear with some machine. 

It is difficult in Case (a) on the scale of our diagram to allow for the second 
harmonics. In Case (@) they become quite insensible and the curves for M@ and M, 
fall together. We may safely say that for y= or <‘01 the values of M, and M, 
are practically identical, and their common value may be expressed by the first 
harmonic 


M=N (1~e, 08 m/ 7): 
41, Now suppose the shaft were to be considered as under statical torsion at 
each instant of time, we should have for the motion of the pulley, rotation ¢: 
Up _ dip 
me a2-—-N+Ks , 
Be being uniform throughout the shaft would = Rees 


Hence we deduce 


where 


Lh 1S SUR eae 
paot— Ns Asin | [ Agee. 


_ aklp oo K ig kl 
Now y= mB and J = Fak — mp? 
if a be here the area of the shaft. 
LN ‘ adi 
Or p=ot—F-+A4 sin Vy + Cf. 
— But o¢=0 and Om, when t=0%*, 


LN st 
; a p= ot—F {1 —cosvy qt 
This gives us 


an -Kt5" =W hi = cosy "7 : 


for all points of shaft. 
Thus we should be taking 


c,=1 instead of 1:00167 (Table V), 
c,=1 instead of °99667 s 
nl=Vy="1 instead of ‘09983 - 


in the case of y=‘01; that is to say, these differences are practically imper- 
ceptible. Or, we may argue that for y= or <‘01 we may treat the shaft as 
taking at each instant the statical torsion due to the couple applied and that 

* Physically of course this state of things could not be set up until ¢=//J, the time the 


wave takes to traverse shaft. This is allowed for in the form given on p. 29. 
8 
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we shall thus obtain a close approximation to the first harmonic term in the 
stress-couple. 
Further we have for the spin 


d IN ee Jt 

deo VvsindyT 
2NS 2 Jt 
a0 AT sin Jy 


This replaces c,='04992 for y=‘01 by //y/2 or ‘05, so that this gives us 
also a very close expression for the variation in the spin. 

Plate XI., Fig. xv1. Case (a) for y= ‘2346 shews how the curves for M, and 
M, begin to separate and at the same time negative values occur for M: as the 
second and third harmonics have a more measurable effect. 

These. results are still more clearly evidenced in Plate VII. Fig. xz. of our 
Art. 22, which is plotted for y=‘5, where we also remark how much more fre- 
quently the maxima begin to follow each other. 


42. In Plate XI., Fig. xv. I have plotted so far as the scale would admit 
the values of M, and M,, and Q—o for y=1. We see here a very considerable 
divergence between the stress-couples at the pulley (M,) and at the rotated 
end (M,) of the shaft. The latter has a very much greater positive maximum,— 
its limiting value being=N (1+ %c,)=2°37N,—and it approaches fairly closely 
to this limit three times before Jt/I=19. Thus for a value of Jt/l between 18 
and 19 we have M,=2:312N. The positive maximum limit of M, on the other 
hand is not closely approached till J¢/] is much greater than the 34 of the 
diagram, not indeed until J¢/l equals about 191, or indeed a better maximum could 
be reached at the double of this. In the shaft 50 feet long dealt with in the 
previous example this is equivalent to saying that the maximum will not be 
reached until 191 x 50/10000 (or even 382 x 50/10000) secs. after the pulley has 
been thrown into gear; this amounts to nearly 1 second (or 2 seconds), within 
which time viscosity has had some influence in reducing the maximum. ‘Thus 
the real maximum-masximorum may lie much nearer the beginning of the oscillations. 

We see further that the negative maximum of MV, is not nearly reached before 
Jt/l=34. M, does not indeed theoretically approach its negative limit before 
Jt/l= about 197, so that again viscosity will have generally had considerable 
influence. JM, on the other hand never becomes negative and approaches fairly 
closely a zero value between Jt//=14 and 15, and Jt/l=22 and 23. The varia- 
tion in Q—wo as here plotted shews only a very slight divergence from the 
first harmonic term. Indeed for all practical purposes for y=0<1 we may 
represent the variation in the spin by the first harmonic. The couples on the 


other hand requiring two, three or even four harmonics to plot them even to 
the above scale. 
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M, should be equal to zero from Jt/l=0 to 1, its slight divergence from this 


value in the drawing is owing to the difficulty of considering a sufficient number 
of harmonics on the scale of the figure. 


We may conclude then that for practical purposes a safe rule would be 
obtained from the following formulae, for y= ors: 
Maximum stress-couple at pulley M,=N {1+ %c,}, 


Maximum stress-couple at rotated end M,=N {1+ Xe} as positive value, 


=N{1—Xc} as negative value. 
RM simu ee: 2IN 
aximum variation in spin at pulley Q—w=+4 e+ radians per second 


_60 JN 


ea ge revolutions per minute, 
T 


where J=velocity of sound in shaft=about 10000 feet per second for a wrought- 
iron shaft, 


N =frictional resistance-couple at pulley, 


K=slide modulus x cross-sectional moment of inertia about axis of shaft, 
when the cross-section is circular*, 


and the values of c,, ¢,, ¢c, are tabulated for various values of y in Tables III, 


IV, V, above. 


* Otherwise the value of KX = slide modulus x vy, where v is given very approximately by 
Equation (iv) of Art. 3. 


Description “of Iuthographed Tables and Plates. 


Table I (a). Values of $(é) and ¢’ (é), Art. 20. Giving effect of throwing shaft-pulley into gear. 

Table I(5). Values of (¢) and ¢’(¢), Art. 28. Giving effect of suddenly applying brakes to wheels 
on an axle, 

Plate I. Fig. 1. Impulsive couple on shaft. Figs. 1. and m1. Impulsive couple of finite duration. 

Plate II. Fig. 1v.. Effect of Repeated loading and unloading, two loads. Figs. vi. and vit. 
Effect of progressive but rapid application of load, Art. 33. 

Plate ITI. Fig. v. Effect of four and six repeated loadings and unloadings on a shaft. 

Plate IV. Fig. vim. Effect of impulsive torsional couples on shafts with fixed terminals. 


Plate V. Effect on shafts when pulleys are thrown suddenly into gear; rotated angles, spins 


al and stresses. Figs. 1x to XI. 
Plate VII. 
eae vt Effect on axles, when wheels are suddenly braked; rotated angles, spins and stresses, 
Plate IX. : 
Figs. xu to XIv. 
Plate X. 


Plate XI. Effect of suddenly throwing a shaft-pulley into gear; variations in spin and maximum 
stresses determined by normal functions. Figs. xv and xvI. 
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INITIAL STRESS IN AXLE WHEN BRAKE /S SUDDENLY APPLIED. 
Vaiues oF Torsionar Couple (KES) FOR FIVE POINTS OF THE 
HALF-AXLE, WHEN Y~/, 2 AND fo. 

ScaLes: Unit of TH - 16 centimetres. 


Torsional Couple. N- 4 cenhmetres. 
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